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Introduction, 

The principal kinetic theories of a gas proceed either on the hypothesis that the 
molecules are rigid elastic spheres, or that they are point centres of forces which vary 
inversely as the fifth power of the distance. Maxwell has worked out the 
consequences of the latter hypothesis in his well-known theory,"^ which is unrivalled 
in its high degree of accuracy and (after some improvements by BoLTZMANNf) in its 
perfection of mathematical form. All the quantities not taken account of in the 
theory (such as the time occupied by molecular encounters, and the effect of collisions 
in which more than two molecules take part) are properly negligible under ordinary 
conditions. The theory has the disadvantage, however, that the underlying 
hypothesis is highly artificial (being chosen chiefly on account of mathematical 
simplifications connected with it, rather than from any physical reasons), and does 
not represent the real facts at all adequately. 

* 'Phil. Trans.,' 1867; 'Scientific Papers,' vol. ii., p. 23. For convenience we shall refer to a gas of the 
type there contemplated as a Maxwellian gas. Of course, its molecules possess no internal energy, 
t ' Vorlesungen iiber Gastheorie,' vol. i. 
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434 MR. S. CHAPMAN ON THE KINETIC THEORY OF A GAS 

The other hypothesis referred to seems to be much more in agreement with fact, 
but its consequences have been worked out less accurately. The method which 
has almost always been used is the one originally devised by Clausius and 
Maxwell ; Maxwell abandoned it later, however, as it had '' led him at times into 
grave error. '^ In spite of its apparent simplicity, numerical errors of large amount 
may undoubtedly creep in in a very subtle way. Hence the theory of a gas whose 
molecules are elastic spheres remains in a rather unsatisfactory state. As a 
'^ descriptive" theory (to use Meyer^s apt term) it has, however, served a useful 
purpose ; the general laws of gaseous phenomena have been developed by its aid in 
an elementary way, which has conduced to a wider diffusion of knowledge of the 
kinetic theorj^ than would have been possible if the sole line of development had been 
by the more mathematical and accurate methods used by Maxwell and Boltzmann. 

In this paper I have applied the latter methods, with an extension of the analysis, 
to the elastic-sphere theory among others. In fact, I have obtained expressions for 
the viscosity, diffusivity, and conductivity of a gas without assuming any properties 
of the molecules save that they are spherically symmetrical. Many known laws are 
thus proved more generally than in any former theories, but the formulae so obtained 
cannot in all cases be put into a really useful form without a knowledge of the nature 
of the molecules. The supplementary calculations required to complete the general 
formulae of Part I. of this paper are carried out in Part II. for three special cases, 
viz., rigid elastic spherical molecules, molecules which are centres of repulsive or 
attractive force varying inversely as the ii}^'' power of the distance, and rigid elastic 
spherical molecules surrounded by fields of attractive force. In Part III. the general 
formulae are completed in these cases and discussed in their relations to the results of 
former theories and of experiment. 



Part I.— General Theory. 

1. Statement of the Problem, 

We shall deal with a gas composed of two kinds of molecules, which are all 
supposed to be spherically symmetrical ; m^ m\ v^ v will denote their masses and the 
number of each kind per unit volume respectively. Similarly the velocity components 
of typical molecules of the two kinds will be denoted by {u^ t\ w), {u\ v\ w^). Let Q 
be any function of the velocity components of a single molecule {e.g., momentum, 

energy). At any point (x, g/, z) let Q be the mean value of Q, so that 

/* (* f 

Q = Qf{u,v^w) du dv dii\ 

where f{u, v, w) is the function which expresses the law of distribution of the 
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velocities among the molecules m.^ Thus the aggregate value of SQ for the 
V dx dy dz molecules m in a volume element dx dy dz at the point considered is 
given by 

2Q = j/Q dx dy dz. 

Let {uq, v^, Wq), {uq\ Vj Wq) be the mean values of (u, v, w)^ {u\ v\ w^). Except in 
considering diffusion, we shall suppose that the mean velocities of each system of 
molecules is the same, so that u^ — uj, Vq = vj, w^ = iVq. We shall Avrite \J = u—u^y 

W = u'-u,\ and so on, so that U = V = W = 0, U^ - V' = W = 0. 

The general equation of transfer of Q, independent of the action of external 
forces ist 



x,y,z 



(1) .^= 2 






here rjrr- denotes the ^'mobile operator" [-t'-^'^o^ — \-VQ-;r- +w^-;r-) of the hydro- 
Dt \dt ax dy ozj 

» 

dynamical equations. The only term of the above equation which needs explanation 

is AQ, which denotes the rate of change of Q due to the molecular encounters ; thus 
the increase in 2Q which is produced in the element dx dy dz by collisions J in time dt 
is dx dy dz dt AQ. The calculation of AQ is the immediate object of our investi- 
gation. 

The motion of the mass centre G of two colliding molecules remains unaltered by 
the encounter, and this point may be taken as the origin of a system of uniformly 
moving axes, relative to which each molecule will describe an orbit in a plane 
through G ; the two orbits will be similar to each other and symmetrical with respect 
to the line of apses. If the molecules move with sufficient velocity to carry them 
beyond the range of each others action, the orbits will each have a pair of 
asymptotes ; the asymptotes of the paths while entering on collision are parallel and 
separated by a distance p (say), and the effect of the collision is to turn the direction 

"^ Thus the number of molecules m in a volume element dx dy dz which possess velocities whose three 
components lie between m and u-\-du^ v and v~{-dvy w and w •+• dw, is 

V f(u,v,w) dudvdwdxdydz; 

this property defines f{u, v, w). Evidently 

/ {% v, w) dudvdw = 1 ; 

/ {liy V, w) is, of course, also a function of a:, y, z and t in general. A similar function /' (i^', v\ w) exists for 
the molecules m', 

t See Jeans' "Dynamical Theory of Gases" (Camb. Univ. Press, 1904), pp. 276-279; or Boltzmann's 
* Yorlesungen liber Gastheorie,' vol. i., § 20. 

I The terms collision and encounter are used indifferently to signify any mutual action of the molecules 
which affects their velocities. 

3 K 2 
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of the relative velocity Vq (say) through an angle 2x in the plane of the orbits (where 
|-7r— X is the angle between an asymptote and the line of apses), the molecules 
travelling away from one another along the second pair of asymptotes with the 
relative velocity unchanged. The angle x is a function of p and Vq which depends 
upon the nature of the molecular forces. 

The velocity components {ui2^ Vi2^ W12) of a molecule m after collision with a 
molecule m! can therefore be written down,"^ from merely geometrical considerations, 
in terms of the original components {u, v, iv)^ {u\ v\ w^), 7?^, m\ x? ^^d e, the latter 
being the angle between the plane of the orbits (which contains Vq and p) and a 
plane containing Vq and parallel with the axis of x. Thus 

(2) u^2 ~ "^-^ 7 [2 {u' —ii) ^^n\—^/{V^^—{u^ — iif] sin 2x cos {e—%vS\, 

m + m 

(3) V,, = v+ -^ [2 {v' - v) 8m\-^/{Y'^{v' ^ vf} sin 2x cos {e^w^U 

m + m 

(4) Wi2 = w-\ — -. [2 (w^—w) sm\—^/{Y^—{w^—^vy} sin 2x cos (e—Ws)]. 

The angle w^ is introduced only for symmetry, as it is zero ; ^V2 and w^ are given by 

(5) (u^^u) {,/-^v) + ^{[Yo'H^'-uy-] [Y,'^{v^-~vy]} CO8W2 = 0, 
and a similar equation in which tv^ — iv replaces v^ — v. Of course, 

(6) Vo^ - (u'^uy+{v'-vy~\-{tv'~-wf. 

This notation having been explained, we return to the consideration of AQ ; we 
divide this into two parts AnQ and A12Q, the former representing the part due to 
collisions of the molecules m among themselves, the latter that due to collisions with 
the molecules m\ Thus 

(7) AQ = AuQ + AisQ. 

Then it is not difficult to provet that 



(* f*-L. CP 



'4- 00 



00 (*27r 



(8) Ai2Q== vv\8y2^)f{u,v,w)f^{%i!,v\%v')dudvdwdu^dv'dw'YQpdpde, 



ao J J J —CO J JO 



where (^igQ = Q12— Q i^ which Q12 is the same function of (t^i2? '^12? ^^12) as is Q of 
{u, V, w) ; thus ^isQ denotes the change in the value of Q for a molecule m, produced 
by a collision with a molecule ml, ^, 

•^ Sea JEANg' *' Dynamical Theory of Gases," pp. 284-288; Boltzmann's 'Vorlesungen iiber Gas- 
t See ihei original papers by Maxwell, or the treatises by Boltzmann and Jea>js, already quoted. 
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The value of A^Q may be obtained from A12Q by making m! = m, u = v, &c., 
throughout. 

The problem before us consists in the evaluation of the integral on the right-hand 
side of the last equation, with certain special forms of Q. Before proceeding with 
this we shall first write down the forms assumed by the equation of transfer (l) in 
these cases. 

2. Special Forms of the Equation of Transfer, 

The equation (l) is more convenient for use in the equivalent form 



V 



DQ _ 8Q " D^o __ 8Q D^?o __ 9Q D^p 

D^ dUa ' Dt dVa ' I)t dWa ' Dt _ 



2 

x,y,z 



i ('UQ)+g{^ (.uv4(«uY)+|(.uw) 



+ AQ. 



In dealing with the problem of heat conduction we shall have to put 

Q = u{u^-\-v^-\- w^) and Q = u^-^-v^ + iv^ ; but in this case it will be sufficient to 
consider the case of a gas at rest, so that u^:=: Vq = w^ = 0, and the left-hand side of 

the above equation becomes simply i^-yt^- ^^ ^^^ right-hand side, in finding 

^^y ^> :t^> we must, of course, not put Uq=^v^ = Wq=^0 till the differentiation has 

du^ dv^ dwo ^ 

been performed. Thus, when Q = u^ + v^ + w^ we have Q = t^o^ + t^o^+t^o^ + U^ + V^ + W^ 
so that TT^ = 7~ = TT^ = 0, and the equation of transfer becomes"^ 

dUa dVa CWr 



/Q ^ VQ ^ M/Q 



Dg _ ^ 3 r 



(9) 3,^^^ 2 ^[,u{U^ + VHW^)o] + Ai,(^^+t;^+t^^), 



where we have written q for |^U^ + V^ + W^; Aii{u^-\-v^-\-w^) vanishes because the 
energy of the molecules, being wholly translational, is unaltered by the encounters of 
the molecules m among themselves. 

In writing out these equations of transfer, it is customary and sufficient to neglect 
the mean values of odd functions of U, V, W in comparison with the mean values of 
even functions, and also to neglect the differences between the mean values of 
corresponding even functions (such as U^, V^, W^) in comparison with these mean 

values themselves ; thus we may write U^=: V^ = W^ = g, and neglect UV, VW, WU 

in comparison with q. Similarly we may neglect UV(U^ + V^ + W^) in comparison 

with U^(U^ + V^ + W^), and the latter we may calculate as if Maxwell's law of 

"^ The suffix on the right-hand side is to indicate that we are considering a gas at rest. Similarly in 
equation (10). 



438 MR. S. CHAPMAN ON THE KINETIC THEORY OF A GAS 

distribution of velocity held exactly. This is because we shall only be considering 
slight disturbances from the uniform steady state. 

Thus from Maxwell's law of distribution^ we find that q = (2hm)~^, 

U^(U^+V^ + W^) = 5gl Putting (^ — u {u^-^-v^ + iv^), and taking u^ ~ v^ ~ Wq — 0, we 
find that t^ = 5g, r-^ = —^ = 0, so that the equation of transfer takes the form 

or (since, as we shall see later, on p. 447, the left-hand side of this equation is of a 
lower order than the terms on the right) 

(10) 5vq^= Au{u'+v'+w% 

We next put Q = u^ and Q = '?^'?; in order to consider the phenomenon of viscosity ; 
in this case, of course, we must take into account the mass velocity of the gas. The 
substitution of these values of Q in the general equation is quite straightforward,! so 
that we set down the results at once : — 

(1.) '*t-^")--«'- 

We now turn to the calculation of At^^ and An {u^ + v^ + w^)o ; the value of Auv will 
not be calculated directly, but obtained from Au^ by transformation of co-ordinates, 
and the value of A^2 {^^ + "^^ -^ "^^ will ultimately be eliminated from our equations, 
and so need not be calculated. 

First of all we shall calculate the values of ^i^Q fr<>na equations (2)-(5) for 
substitution in (8). Since, however, several terms of ^isQ disappear on integration 

with respect to e, it will be most convenient if we immediately calculate S12Q d^. 



r2n 

3. Values of ^iJ^de 

Jo 



In calculating S^J^ we shall find it convenient, partly for immediate brevity, but 
much more for the sake of a subsequent transformation of the variables of integration 
in (8), to write 

(13) (m-fmO (U'-U) = X, {m-^m') (V^-- V) = Y, (m-fmO (W- W) = Z, 

(14) mU-f-m'U' = X^, mV-f-m'V^ = Y^, mW + mW = Z^. 

* In which fill, V, w) - /^y^' e~-hra%-v^ • 
t See Jeans' treatise, §§ 336, 338. 
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Thus 

and 

(15) (m + mO (mSU^+m'SU^^) = Vi^+mm%^ 

writing Y,' = X/ + Y/ + Z,'. Also, we have 

(16) (m + mOU = Xi-FX, (m+mO V = Y,^/c% (m+mO W = Z^-y^'Z, 

(17) {m+m')W = X, + kX, {m+m')Y^ = Y, + kY, {m+m')W - Z. + kZ, 

where A; = : , k = 



Since the mass velocities of the two systems of molecules are the same, we have 

{m-\-m!){u^—u) = (m + m^)(U^— U) = X, 

and two similar equations. 

Considering Q = u^^ it is easily seen that 



'2w 





C2ir 




{Ui2—U^) de 



m! 



= / _^ A2 [(^'-'?^) (m^ + mV) Stt sin^ x+^' {Yo'-3 (^'-'?^)^} ir sin^ 2x], 

where we have expressed the result in terms of sin^ x ^^^ sin^ 2x ; on making the 
above transformations we obtain 



'^"' m! f-xr -XT . / . n ■xr■^ . • o . m^^ 



(18) S,2u'de = j-^^^^-j^iX,X-h{m^ 

We next consider Si2u{u^+v^-i'W%; we write t^ig = '^ + <^— <^^ cos (e— t^i), where 
a = 2F (u^—u) sin^x? ^^ = k'\/YQ—{u'—uf sin 2x, and similarly for i^ig and Wi2' Since 

'^i2^+'^i2^+'2^i2^ = 2 ('U + a)^-2S ('U + a) a^ cos {e-w^)-\-i:a^^ cos^ (e— 'i^i), 
we have 

— u^2 h^ii + '^12^ + ^^12^) de = ('2^ + a) 2 {u + a)^ + 1 ('u + a) Sa^^ + a^S (t^ + a) d cos 'i^;!. 

-^TT Jo 

It is easy to show that 

a'S (w + a) a' cos w^ = , ^'\.^ {u iY'+Z')~vXY-wXZ} sln^ 2y, 

and from the last three equations it readily follows that 



*^ ml 



(19) .Jl3«(^*Hv^+^<;^)oC^e= 7-^^,(2Xi2XXi + XV/+m'^XV/)4xsin^v 

Jo (m + w; 

2m'^ 

— trrcSXSXXi — XiSX I TT sm 2y. 

(m+mO^ 11/ A 
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Integration with Respect to p. 

The integrand of equation (8) contains p explicitly in the term p^ and implicitly 
in X (which is a function of p and Vq) ; the latter occurs only in the form sin^ x ^^^ 
sin^ 2x. Hence to integrate with respect to p we need only know the values of the 
two integrals^ 

47rVo sin^ X • P ^JP? '^'^o sin^ 2x . p dp, 

Jo Jo 

which are functions of Vq only. It is sufficient, at present, to denote them by the 
symbols O'(Vo) and Q^^(Vo), leaving their further consideration till later. As the 
functional relation between x? "^o? ^^d jp, and consequently also the values of 
the above integrals, depends on the law of molecular interaction, which will differ 
according as the collision in question takes place between two molecules m, two 
molecules m/, or one of each kind, we distinguish between the three cases by adding 
the suffixes 11, 22, 12 respectively to Q^ {Y^) and Q^^(Vo). 
Thus 

(20) A,,w^ = r~^^ \\\[\ [2 {XX, + (m + mOt^oX} Q'^Y,) 

\fjfi -T ni/ J J J J J J J 

+ F ( Y^ + Z^-- 2X^) Q^\, ( Vo)] f{u, V, w) f {u\ v\ w^) du dv div du' dv' div\ 

(21) A,;^ {u^-^v^+w% = ^j;^^,^, WWW [(2X,2XX, + XV/+m^^XVo^) Q\, (Vo) 

In the case of a Maxwellian gasf the expressions (Vq) are independent of Yq, so 
that the integrals just written can be evaluated in terms of mean values of functions 
^f U, Y, W without any knowledge of the functions f{u, v, w), f {u\ v\ tv^). In 
general, however, this is not possible ; we require some knowledge of these functions, 
which express the law of distribution of velocity, in order to make further progress. 

4. The Laiv of Distribution of Velocity. 

It is well known that in a gas which has had time to attain a uniform state the 
functions /, f^ have the respective forms 

hm\ ^-hm^jj^- ( hm \ ^-^hm'%v"' 

IT J \ TT 

* The factors 4:7rVo and ttYq are added merely for convenience. 

t Also in the case of a gas whose molecules are point centres of force which attract one another 
according to the fifth-power law. This case, among others, is considered in Part II., and probably is 
nearer to actual fact than Maxwell's case. — Note added October, 1911. 
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TE^'X' 



which we shall denote by /o? /V When there are slight inequalities of temperature 
and mass motion in the gas, we shall suppose that 

where F and F^ are expansible in the form of power series (without the constant 
terms) in the variables (U, V, W) and (U^, Y^, W^) respectively. The coefficients 
will be small quantities which are functions of the velocity and temperature, and 
their derivatives, at the point considered. Thus F and F^ represent the small 
disturbances from the normal law of distribution, caused by the slight lack of 
uniformity in the gas. Maxwell and Boltzmann'^' considered that the terms of 
the first three degrees are sufficient for the adequate representation of the disturbed 
state, and I shall follow them in this assumption. Thus we write 

(22) F (U, V, W) - {2hmy" {a,'U+a,Y+a,'W) + 2hm (i^anUH 22ai,UV) 

+ {2hmfl' (-i-2aniU^ + i^a,,,VY' + a,,,VYW) 

where the factors {2hrrif^^\ {2hm), {2hmf'^ are added merely for convenience in the 
integration. We shall have a similar equation for F', in which m, ^i, ag, ... are 
replaced by m\ a\, a\, .... Since, by definition, 



f{u^ V, w) du dv dw = 1, 



J %> J 



f {u\ v\ w^) du^ dv' dw' ~ 1 



we have 



or 



l + aii + a22 + (^: 



33 



1 + ^11 + ^22+^33 ~ Ij 



(23) 



an + %2 + <^33 = 0, a\^ + a! 22 + o!zz = 0. 



* By Maxwell in his memoir "On Stresses in Earefied Gases/' *Phil. Trans.,' 1879, or * Scientific 
Papers,' vol. ii., p. 681 ; by Boltzmann in his ' Yorlesungen iiber Gastheorie,' vol. i., p. 185. In each case 
the assumption was made in connection with a Maxwellian gas, but there are good grounds for believing 
that it is equally valid in general. As neither Maxwell nor Boltzmann considered it necessary to give 
any justification for their procedure, I deliberately followed the same course, more especially as the 
attempt to make the step perfectly rigorous would have necessitated the introduction of much mathe- 
matical analysis which would be out of place here. 

I should also mention that Enskog ('Phys. Zeitschrift,' xii., 58, January, 1911) has made an attempt to 
determine directly the form of the function F iii, v, w), applying methods of integration, similar in many 
ways to those used in this paper to evaluate AQ, to an equation arrived at by Boltzmann (' Vorlesungen,' 
vol. i., p. 114). From the expression for F (u, % w) thus obtained Enskog deduces values of the coefficients 
of viscosity and thermal conduction for a simple gas. 

I am indebted to Prof. Larmoe for the reference to Enskog's work, of which I was unaware till after 
this paper had been communicated to the Royal Bodety.—Note.. added Odoher, 1911. Some of these 
statements are modified by the last note on p. 483. 
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Similarly, since U = Y = W = 0, U' = V = W = by definition, we have 

(24) 2(*i + aiii + aj22 + ai33 = 0, 2a'i + a'ui + a'i22+a'i33 = 0, 

and foTir similar equations.* 

We require the mean values of the following functions of U, V, W : — 



U^ = {2hm)~^ (1 + fflii), UV = (2Am)-i a,^, W = {2hm)-'l' a,j„ UV^ = {2hm)^''' a,,, 



1225 



similarly for the second system of molecules. 

As usual, being the absolute temperature, we have {2h)~~^ = E0. The various 
components of partial pressure due to the first gas, p^^., p^^, &c., are given by 



»U' = ^ (1 + an), p,n = /)UV 



V 



p,, = pv- ""2h^'''' ^'' ^ ^ ^2h'' 



since p ^ vm\ the mean hydrostatic pressure p is given by 



similarly for the second set of molecules. 

In substituting j5^^ in the equations (20), (21), we shall write 



TT 



(by equation 15) since FF^ being the product of two first order small quantities, is 
negligible. 

■^ It should be noted that the above expression for F is of the lowest degree consistent with the 
satisfaction of the requirements. The function F must provide for small changes in the mean values of 
even functions of U, V, W, and also of odd functions, both these changes being of the same order. The 
terms of the second degree do this for the even functions, an, ^22, %3 being of the first order; it might afc 
first sight be thought that the terms of the first degree would provide for the odd functions of U, Y, W, 
but this is not so on account of the conditions U = V = W = 0. Hence the terms of the third degree 
must be present, and their coefhcients must be of the same order as %, an, &g,— Note added October, 1911. 
§ee the note on p. 483, 
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5. Reduction of the Integrals. 

The evaluation of the expressions (20) and (21) is facilitated by changing the 
variables (U, V, W), (U^, V^, W^) to (X, Y, Z), (X^, Yi, Zi) in accordance with equations 
(16) and (17). The Jacobian of transformation is easily found to be (m + m^)"^, and 
the limits remain as before, viz. , — oo and + oo . The two integrals are now of the 
form 



(* r 



A 



J J 



^-«2x^-5sx.^ i^ ( Vo) V>^ (X, Y, Z, X,,Y,,Z,){l-{-¥+¥)dXdYdZdX,dY,dZ, 



where \fr is an integral polynomial in the given variables ; F and F^ are also 
polynomials in these variables. Evidently only those terms in the product of ^r and 
(l + F + F^) which are of even degree in each of the six variables separately will give 
any result upon integration. It is easy, though tedious, to pick out these terms ; 
evidently if \fr is of odd degree in the six variables combined, only those terms 
of F and F^ whose a- coefficients have an odd number of suffixes will need to be 
considered ; similarly, if i/^ is of even degree, we need only consider terms in F which 
have a^, a^s, ... for coefficients. Having picked out these terms we are left with a 
number of integrals of the form 



/* I* 



I* f* 



%/ *j */ */ 



i^-a^x^-,tx? Q (Y^) X^T^-Z^%^^'Y/-Z/^' dX dY dZ dX, dY, dZ,. 



The integration with respect to Xi, Yi, Z^ can be carried out immediately — ^most 
conveniently by changing to polar co-ordinates. We do the same also in the case of 
the variables X, Y, Z changing to the variables Vq, 6, (p. The integration with 
respect to the latter two variables is simple, and there remains an integral of 
the form 

where a = — ; — ^, ^hm, or ^hm! according as Q(Vo) or A has the suffix 12, 11, or 

22 respectively. 

We shall consider this more particularly later ; at present it is sufficient to denote 
it by a symbol. As we shall only need to consider the three cases ^^ = 2, ti = 3, ti = 4, 
it is, perhaps, most convenient to denote it by P, E, S respectively, instead of 
adopting a more general notation. To distinguish between the different cases arising 
from the various functions Q (Vo) occurring in the integrand, we add the same dashes 
and suffixes to P, R, S. Thus we have integrals P^g, Pn, P^n, and so on, 
corresponding to the cases when in the above integral ti = 2 and Q (Vq) has the 
special forms Q^^y ^ n, ^^i, and so on : similarly for R and S. 

Though the execution of the processes indicated is rather lengthy, it is quite simple 
and straightforward, so that, without entering into the details of the calculation, we 

3 L 2 
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shall quote the results forthwith. After simplification by means of the equations 
(23) and (24)5 we obtain the formula 



(25) A^,v/ 









2 I ^ \ '\^r..^Jr.. I ^1 



a\ 



JlJ 



^^^3/2 ^/8/2 



/ 2 

12 "~ \5 



o hmm^ 



m-\-ra 



f '^ v. 



+ l(%i-OF 



_U %. ^™^^-'^ /^ll ! /-^lM ][>^ 



2 
12 ' 5 



-^11 . ^ 1 1 



7/ 



Qn-f-m \in 



m 



(26) Ai#K'^^ + '^'^ + '^^^)o = 47r"^^V 



mm^^ / hmm^ Y^^ f2Y'^ 



(m -l~ m^)"*^ \m + quw \h/ 



_8 

1"5 



a-t 



a 



hmm!'^ 



m-\-m!\m^^^^ m^-^^^, 



9. ^ n h> J-^ 12 



I r}%w}! I 1 



6 



/., + ^ /i::i;U>^ ^'^ 



m 



\wi' 



3/2 



m 



/3/2 



{m}^^ai-~'7n!^^^a 



1/ ( JL 12 



where for convenience we have written k^ and /{.^ for the quantities given by 



(27) 

By putting 
following : — 

(28) 
(29) 



7 _ 2 h'in'}n^ Wj 
m + 'in r 



2 

! > 

12 



h 



2 hmm! \^ S^^. 



m + m 



12 



m 



m^5 1/ ^ 1^^, a = a! in the above equations we get also the 



Aiit^^ 



%rrb\^^^ hm 



-47r.M^) ^C^.R^^ 



(TT 



20 



11-^^ IV 



^{li{u^ + V^ + %V^)Q = 4.Trv' 






i^ a, {2hmy" W\,. 



By transformation of co-ordinates we obtain the following equations from (25) 
and (28) : — • 



(30) A 



-~47r pv 



mm 



f2 



Jimm' Y^'^ 



{m-\-'iv!Y\m-\-m^, 



hmmJ^ fa 



|(m+mO(ai.-agF..+ 1 '£^, (^^ + ^]^"v- 



+ 1^ {2}if'^mm^ \v 



a. 



a\ 



^^^3/2 



m. 



/8/2 



+ % 



a 



Vn 



a 



! \ 



n^^l2 



m 



/3/2 



F -2 



7n-\--m' \7n 
hmm! 



9l 

m,' 



12 5 



m^ + w. 



/ J^ 12 



(31)^^^ 



A.MV 



^11 



4zlTV 



limY^^ hm 



2 



TT 



20 



^12-^ 11' 



^ All the formulge (25)-(31) were in the present form before the paper was revised, with the exception 
of (26). Their calculation was given almost in full, and was performed by essentially the same methods 
as those explained above. The work was made rather more lengthy than necessary by the use of an 
unsymmetrical transformation of the integrals in place of that given in equations (13)~-(17), to which I was 
led while endeavouring to simplify the calculation of Av^i{tt^ + v'^ + ^2)0. The calculation of Au{to^ + ^24. ^2^ 
was also complicated by the fact that the gas was not assumed to be at rest; but on revision, despite these 
simplifications, I decided to omit the routine calculations altogether. 

Equation (26) was not given in the paper as originally written because it is connected with the 
conductivity of mixed gases, and I was then unaware of any experiments on the subject which would 
make so tedious a calculation worth the while. I. have since found stich experimental data, and have 
therefore worked out the formula (26), the results of which (as will be seen in Part III.) show very 
satisfactory agreement with the observations. — N'ofe added October^ 1911. 
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6. Calculation of Ayju. 

Finally we turn to consider A^a'?^^, the calculation of which is rather different from 
that of the foregoing quantities because Uq and u^^, v^ and v^q, w^^ and iv^^ are now" 
different. We require Ily^u in order to determine the coefficient of diffusion. 

It is evident from equation (2) that 

§12U de = 47r : iu^ — u)^ 

m + m 



so that 

A / ^' 

l\l2U = VV . 

m -h m J 



J J J 



0^2 ( Vo) (u^—u) f{u, V, %v) f {%\! , v\ tv^) du dv dw du^ dv^ div\ 



In the present case we may neglect the small deviation from Maxwell's law of 
distribution of velocities, and thus avoid making any assumption as to the form of F. 
This is so because in this case (unlike all the preceding cases considered) the unit 
term of (l+F + F^) in the expression for ff^ leads to a term in the final result, so that 
the terms which would result from the inclusion of F and F^, being multiplied by the 
small coefficients a and a\ are negligible to our order of accuracy. The term arising 
from the unit in (l + F + F^) is itself a multiple ot u^q--Uq^ which we assume to be so 
small that its squares and higher powers are negligible; ii^q—Uq being thus a small 
quantity of the first order, the terms arising from F and F^ are of the second and 
higher orders of small quantities. 

Thus we now have 



2 

IT 



p~h [im-bm') %it^-\-2^ (muo+m'u'o+x) u—m't (w'o+a;)^-wSwo'^] 



where we have written {x, y^ z) for {u—u^ v—v\ iv—w^)^ so that Y^ = x^-{-y^-\-z^. 
Changing the variables in the above integral to {x^ y^ z) in place of u\ v' ^ %v\ and 
recalling" that 



»Q0 



00 






on integration with respect to u^ v, iv we obtain 

Ajgt^ = —VV -, i —7 }\ r 12 (VQJxe '''^+''' dx dy dz 

m-hm [7r{m-{-m)} jjj 

where we have written (Vo^— 2yoVo cos X + /)^) for 2 {x-{-u^Q—UQy ; thus 
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Since, as has already been mentioned, p^ is a small quantity of the second order 
the term e ^^+^*' may be neglected, since it is equal to unity to our order of approxi- 



mation. Transforming the variables x^ y^ z in the last integral to polar co-ordinates, 
we have 



Aiot^ = — j/j/ 



m 



/ 



hmm! 1 ^^^ 






ra-\-rrh^ l7r(m + m^). 



hmm' 



Vo^0^12(Vo)e «^+«^' 



7V0' 



2 



y 1. (hmm! 



7 2pVo cos X ) > (iVo sin cos Q dQ dip 



n! \m-\-m[ ! J 



by putting x — Vq cos 0, &c. It must be remembered that X is a function of and ^. 
since 



./ ., .y 



cos X = -^— — - cos 0+ -^^ sm cos ^ -{ — ^^ ~ sm sm ^. 

p p p 

All the terms of the exponential series occurring in the integrand of the last 
integral are negligible (on account of the factors f?^ yo'^, and so on) except the unit 
term and the term of the first degree. The unit term leads to a null result on 
integration with respect to and (p. The first degree term alone contributes to the 
final result, which (in our previous notation) is easily seen to be 

(32) Ai,^i = |x-n/-^(^^,'f (<~«o)P'i. 

^ m-\-m \m-\-m J- 

If we put m = ml, Uq = u^q, j^ = / we find that A^u = 0, as is otherwise evident, 
since the momentum of a system of molecules is unaltered by their mutual encounters. 

Having now calculated all the values of AQ which we require, we proceed to 
substitute them in the various special forms of the equation of transfer, in order to 
obtain expressions for the coefiicients of viscosity, diffusion, and conduction in simple 
and mixed gases. 

7. The Coefficient of Viscosity of a Simple Gas, 

First considering a gas composed of molecules of one kind only, we substitute for 
Att^ (which in this case equals A^u^) and Auv from equations (27) and (30) in the 
special equations of transfer (ll) and (12) respectively. We get 

'"^W 2F^^^^^^-"~^na^"'¥, 

Remembering that a^ = — {Pxx—p)^ ^12 = — Pxy (see p. 442), and comparing these 

V V 
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equations with those of a viscous fluid whose coefiicient of viscosity is /x, viz.. 



Pxx-P = M 



_2 
3 



dx dy dz / dx. 



^^^ ^ \dx dy 

we obtain the following expression for fx : — • 

(33) 5(7 (2.Y 1 - 5 ^^^r ^ 



8. J^^e Coefficient of Thermal Conductivity of a Simple Gas. 
We next substitute for Au {u^+v^+w% in equation (lO), and so obtain 






on substituting for aj from the equation on p. 442. This equation is simplified if we 
recall the value of /u. from equation (32). Thus* 



3^=-ltU(U^+V^+W^,). 



If we now substitute for U(U^ + V^+W^) in the equation (9), we get 



Dt x,y,z'dx\ ^ dxj' 



If we compare this with the equation of conduction of heat in a gas at rest, whose 
thermal conductivity is 3- — viz., with 



5^ = 2 —1^9-—' 

J^t x,y,zOX\ CX/ 



(where C^ is the specific heat at constant volume) — and remember that q is directly- 
proportional to the absolute temperature 0, we find that 

(34) ^ = 1^0,. 

^ We can now see why the term v — U (U^ + V^2 ^ ^2) ^^ the left-hand side of the equation preceding 



equation (10) is negligible. The ratio of its coefficient v to the coefficient of U (U2 + V^^ W^) on the right 
hand (viz., in ^u {ii? -{• v'^ -\- vF)q) is now seen to be /x/p, which is an exceedingly small quantity in all gases 
Under normal conditions. 
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9. The Coefficient of Diffusion, 

The general equation of transfer (l) is true for any system of molecules in a gas, 
whether in the presence of other systems of molecules or not, provided that no 
diffusion is taking place (this restriction arises in the elimination of the external 
forces,'^ where it is assumed that Hu =■ 0, which is true, of course, only when the 
restriction just mentioned is satisfied). The form of the equation of transfer which 
is applicable to the more general case now under consideration, where diffusion is 
taking place, isf 






dt 



X, y, z 



cx ex m cUo_ 



+ AQ. 



where X is the component of external force, If we put Q~u, and suppose that there 
are no external forces, and that the motion has assumed a steady state, the equation 
becomes 

ax 

where we have neglected products of u^^ Vq, iVq as being small quantities of the second 

order ; similarly we have neglected such terms as UV, and have given i/mUQ its 
proper value p. 

The equation just obtained applies to the first system of molecules, p being the 
partial pressure due to this system ; there is a similar equation 

7^ = m^A2iU 
dx 

for the second system. Since the temperature is supposed to remain constant, we 

have 

dp _ _3_. /_j^\ _ J_ _8i; Sp'' __ J_ ^ 

cx cx\2hj 2hdx^ dx 2hdx' 

Now, from equation (31), it follows by symmetry that 

m^A2iU ~ —m-Ay^u^ 
so that 

;^ ——■;;— -r -§- i/^ ■ ■ ; ; / TT All [U q-~U^} iT y^' 

CX CX ni~\-nh \in + m/ 

The total flow of molecules of the first kind per unit area per unit time is clearly 

vu^y a.nd also (by the definition of the coefficient of diffusion D^^ is equal to — D12 — ^. 

cx 

^ See Jeans' treatise, p. 279, 
t Ibid., p. 278 
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Hence 

cx dx ex 

by which the preceding equation may be reduced to 

dv 



RoD 



12 -■-'12 



dx 



dx m+m xm+m / 

This gives us* as the expression for D12 

By putting m = m\ v — v\ we get the following expression for the coefficient of 
diffusion D^ of a gas into itself : — 

(36) Dh = -/^Y''(M''^'-- 



V \2i X 11 



11 



10. The Coefficient of Viscosity of a Compound Gas, 

As in dealing with the case of diffusion, we must now use an equation of transfer 
for each system of molecules. Writing 

dx ^ \dx dy dz / ■ ^ dx dy 

where (since the various systems of molecules are not supposed to be diffusing 
through one another, so that Uq = %i\, &c.) E and F are the same for both systems of 
molecules, by equations (ll) and (12) we have 

In the expressions already obtained for Ajg^^, l^i2^v (viz., (25) and (29)) there occur 
terms containing the coefficients a^, a\. As we have already seen in discussing the 
conduction of heat, these coefficients depend upon the existence of variations of 
temperatures in the gas. We shall here suppose that the gas is at a uniform 
temperature throughout, so that the said terms will disappear. 

* Since this paper was written I have found that the expression (35) had already been obtained by 
L ANGEVIN (* Ann. de Chimie et de Physique,' (8), v., 245, 1905), who applied it to the motion of electrons 
in an electric field. The present proof is shorter than Langevin's. Enskog (* Phys. Zeitschrift,' xii., 
533, July, 1911) has also published a simplified proof on lines not unlike those of the above proof. I arn 
indebted to Prof. Larmor for the reference to Enskog's paper (which appeared while this paper was in 
the hands of the Eoyal Society), where I found the further reference to Langevin's theory.- — Note added, 
October, 1911. 
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Thus we have 



mm 



^^ I hmm! Y'^ 



{m + m')'^ ytn + m' 



2T>^ I 2 ^^^ "R// I 



;J: TO + 



-.R' 



m + m 



a 



11 



1 



2P^ 2_ 

3-^ 12 5 



" r 7 -tt 12 f <^ 11 

m + m J 



If we write 
(37) 



k = ^^ 



5 



hmm! 

m + nh 



.15// W2_p/ 

/ J-^ 12 // 3-*^ 12) 



and substitute for P^jg i^ terms of the coefficient of diffusion Dig, the equation for 
Ai2U^ is simplified to 

1 



''^12 



1 ^^ 

2 



m 



V 4- / K^mmf {m-\-m') Dj, 



1 +/t 



\j 



m 
m 



^n— i — 



ai 



il-k)a\^> 



Again, by means of the expression obtained for /ul, the equation (27) for Autt' 
reduces to 

11 



Ante^ = — 



v^a^ 



'11' 



4h^mjn * 



We get similar equations for Ante^^ and Aigt^^^, writing jul^ for the coefficient of 
viscosity of a gas composed only of the second kind of molecides. Hence^ 
remembering that AQ = AnQ + A12Q, we have 



E = Aaii + Ba^ii, E = Can + Da' 



where 



XX 



B = 



C 



D = - 



V 



1 + k 



/\ 



11 5 



V 



V + V h{m + m!) D12 
1 1 



m 



2h 



/A 



V 



V + V h (m + m^) D12 



p 



p + V h{m^-Vm!) D^g 
' 1 1 



(l-Zc), 



p 



1 4-/ — 
p + p^ h{m + m^) D70V m^ 



V 



JjfhfJ> 



We next substitute for an and a\i their values — _p^,,, tP^xx respectively in the 

p p 

above equations, and solve so as to obtain P:,:,+p^xx (the total normal pressure parallel 

to O^) in terms of E. After a little reduction we find that 



jPxx "T"jr XX 



A'+B'+C'+D' 



A'D'-B'C 



/n/ 



E, 



where 
A' 






P \ Wj IK 
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In these last expressions w and w^ are the specific gravities of the two component 
gases referred to a standard gas whose density at the pressure and temperature of 
the compound gas here considered is p^ ; and G is given by 

(38) G= ^ ^ ^ 



Since the partial pressures of the component gases are directly proportional to v, v 



f 



for jt> = -^ , ^' = ~7 ) , we may write 



Similarly we may show that Pxy-^P^xy is the same multiple of F. Hence, recalling 
the denotation of E and F, and comparing the said equations with the equations of 
pressure in a medium whose coefficient of viscosity is /x^g? ^s in § 7, we obtain the 
following equation : — 



W ^ ,2 



^'' 1 /1 . 7 ^^\ 2 , {{W±wj .^ , 11 ,1 /' , ,tv'\ r, 

M \ w J L WW (jTyu/x J /x \ ^(; / 



11. 2%e Coefficient of Thermal Conductivity of a Compound Gas,^ 

From equations (26) and (29) we see that /^u{uF'-\-v^-\-%v^)^ and Lu' {u^^-{-v^^-\-w^^^)^i 
can be expressed as in the following equations : — 

i^u{u'^v'^-w% = -2v {A(2Am)-^/^ai + B(2AmO-'^'a'i}, 
Lu^{u''+v''+w% = -2/{C(2Am)-^/^ai + D(2AmO--'^'a\}, 

where A, B, C, D are written for convenience in place of some rather long expressions 
which can easily be written down from the equations cited. 

Remembering the values of a^ and a\ as found in § 4, and substituting from 
equation (lO) for I^u{u^-{-v^'\'W^)q^ the last equations may be written 

5.g|2= . (A U(UHV^+W^)o +B tF(U^^ + V'^+ W'^)o), 



5.V ^ = / (C U (U^+ V^+ W^)o +D U^ (U^2^ Y/2^^/2)^)^ 

* Added October^ 1911.—- See the footnote to p. 444. 

3 M 2 
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Since q = 1/2 Am, so that mg = m^q\ on solving for U (U^ + V^ + W^)o we obtain 

(AD-BC)U(U^+V^+W^)o = 5(^D-^b)-^, 

\m m I ox 

(AD-BC) U^(U^^+V^^+W^^)o = bl^K-^G)^. 

\m m ■ I ex 

If we multiply equation (9) by m, and add the corresponding equation for the 
second system of molecules, we get 

* 

since the sum of the remaining terms 

mAi2 {u^ + v^+w% + m^Ais {u^^ + v^^ + iv^% 

vanishes. This follows from the principle of conservation of energy, for the last 
expression represents twice the rate of change of the combined energy of the two 
systems of molecules due to their mutual collisions, which is evidently zero. 

We substitute the values already obtained for U (U^ + V^ + W^)o and 



U^(U^^ + V^^ + W^^)o in the last equation, and compare the resulting equation with 
the equation of conduction of heat in a medium at rest. As in the case of a simple 
gas (§8) we obtain the following expression for the thermal conductivity S-jg : — 



^- — 5 



m m m m 



,+/ AD-BG 



In this equation we must substitute the proper values for A, B, C, D. I shall not 
enter into the details of the calculation, which is rather long and complicated, but 
will simply quote the result in the simplest form. It must first be mentioned that 
(C^)i2? the specific heat of the compound gas at constant volume, is connected with 
the same constants for the component gases, viz., C„ and C^^, by the following 
relation : — - 

(39a) (a),. = ^^^-+^yp^- 

pw -\-pw 
using the notation of § 10. 

The formula finally obtained is 

^ ^ E — J9^ + Fi PJ9 -+ G —,%> 

fj, fJ. 
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where E, F, Fj, G are given by the following equations, in which the quantities 
h, ki, h^ have the values assigned to them in equations (27) and (37) : — 



r 



(41) \ 






G 



P = 



%lc 



w . w 






w w 



2 {w-^Wy J. (w 

WW • \/>f 



+^)-j-¥^+*^'-^\i-^0-Hi + 5^.) 

yU / L ' WW : 



V, 



(w + tt; ) /)o-Ui2 I , WW , , \w wV J m/ 



D 



12' 



This completes the first part of the paper. We proceed now to the evaluation of 
the quantities P'12, Ei'ig? ^j ^n? ^^^ ^2 '^^ some special cases. 



Part II. — On Certain Special Forms of Molecitlar Interaction. 

12. The general expressions which we have obtained in §| 7-11 contain four 
integrals, P^g? ^^2? ^'^2? 8^2? given by the equations 



p/ 



(42) 



•00 



hmmf 



12 



Vo^Q'i, (Vo) e-^' '" dVo, 



E" 



12 



•CX3 







hmm' 



^Vo^ 



Yi^"n{^,)e~^^-'° dY„ 



< 



R',2 = r Vo«0',, (Vo) e 

Jo 



hmm/ 



m+m 



-.Yo' 



dY„ 



S' 



•00 



12 



Vo«ii\, ( Vo) e 



hmm' 



m+m' 



}W 



dV.. 



0' 



The remaining integrals, P'n and W\i can be obtained at once from P'12 ^^.d W\2 
by changing m to m^ and 0^2 to Q^, 

We recall the values of Q\2 (^0) ^3:id ^'^2 C^o) from the latter part of § 3. They are 



(43) 



•00 



^^2 (^^o) = 47rVo sin^ X • P <^P> ^'^2 (^o) = '^^0 sin^ 2x*pdp 

Jo ~ Jo 



'00 



where 2x is the angle through which the direction of the relative velocity Vq is 
turned by the mutual action of the two molecules 7n and mV while p is the distance 
between the ^.symptotes of their orbits relative to their mass centre (§ l). 

As the above expressions cannot, however, be integrated except in certain special 
cases (which fortunately include the most interesting cases) we consider these in 
order. 

13. Rigid Elastic Spheres. 

The simplest case is that of molecules which are rigid and perfectly elastic spheres, 
of radii a- and o-V Clearly it p exceeds o^+o-' the spheres will not affect one another's 

* ' 

motion, so tfeat x = 0- Whenp is less than o-+a-^ it is evident that 

P 
smx = ~~f-r 

c-ro' 
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Hence 

^0 \ar-\-ar ) 

and 

r(T-\-<T' 

^'^12 (Vo) = ttVo ((t+o-^^ sin^ 2x sin x cos x^X- s'ttVo (a- + crO^ 

Jo 

Consequently we find that 

(44) V\, = ^ (<T+<7'? , Yi e"^'^ " dV, = TT (^+crO^ (W^ , 

Jo Xhrnm / 

A4 



\ hmm I 



y\5 



S'i2 = 127r (o' + o-^)^ ( -T- J ) • 

\ hmm I 
From these we obtain the equations 

(45) P^H = 32w^ f^Y, R\ - 64' *^ ^ ^ '' 



2 
TTCr 



so that in this case the constants ^, /^i, ^2 (defined by equations (27) and (37)) are 
numerical constants ; in general^ they are functions of the temperature (^'.e., of li). 

14. Molecules which are Point Centres of Force, 

We next consider the hypothesis that the molecules are geometrical points 
endowed with inertia and repelling or attracting one another with forces which are 
functions only of the distance between their centres. 

Let ^12 (^) be the mutual potential energy of two molecules m, m! at distance r, and 

let us write ^^12 (^) i^^ place of j- ^12 i"^)- The first two integrals of the equations 

of motion of the second molecule in terms of co-ordinates r, 6 with the first molecule 

as origin are as usual 

« • 

r^0 = A, ^{P+rW) = ^^12 +B 
where A and B are constants. Eliminating the time from these equations we get 



*{(IJ--}"=*'"«-B - (IJ=S 
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Taking the direction of one asymptote for initial line this has the integral 



e 



GO 



2B 






1/2 



dr. 



Since Vo is the relative velocity before collision, and j) is the perpendicular from the 
origin on to the asymptote, we have 



so that 



A=i9Vo, B = iVo^-/,,(a3), 







CO 



r 



4 



Lf 



,^+ 



2r' 



A-fT <^ 



{<p\2{r)-i>\2{^)} 



-1/2 



dr. 



The apse is given by the vanishing of the expression in square brackets. If ro is 
the positive root, we have 



X 



•ro 



00 



r 



.4 



■P 



2 



'^'2_j_ 



2r' 



p'Y.. 



2 



^'i2{r)-<p'i2{^)} 







-1/2 



dr 



since x is the angle between the asymptote and the apsidal radius vector. 

For some purposes it is more convenient to transform the variable of the last 
integral to j;, where fj = p/r. We thus get 



X 



'% 







2 L' fp' 

iiMimiimriiiiiwint^Hp- <1 #T\ I ■^= ^ —»w 



l-'?'+^i<i^'l2(^)-<^'l2(«=) 



-1/2 



C?;?, 



^0 being the least positive root of the expression in square brackets. 

No further progress can be made without knowledge of the form of ^'12 or ^12- The 
simplest and most natural form to consider is that corresponding to the case in which 
the molecules attract or repel one another with a force varying inversely as the 
fii^^ power of the distance. In this case 



^12 (r) 



K 



CO 



/v»?il2 



- mm^ dr 



mm 



K 



^^2('^) = —{m+m^) 






(%2— 1)^-"^'"^' 



where K is a constant depending on the nature of the molecules. 
Hence 



X 







^ 



{t^-^l)VoV" 



•1/2 



d 



n 







1 ^ ' "^ 



n~l 



n 



.a, 



•1/2 



d 



n 



where"^ 



a =J9 



n 






2Ki2{m-^m^). 



7132—1 —y. — J— , 

V «12-1* 







Thus X is now a function of a only. 

* When the forces are attractire the sign of a^""^ is reversed. In the text the case considered is that 
ol repulsion. 
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We now transform the variable in the integrals Q\ Q" from p to a ho that 



where 



Thus we have 



pdp = K^i2^o ''^'"^ ^ da 



-•-^12 — 1 T 

I %2"~1 



?ll2--l 



0',,(Vo) = K',,Vo^-'^ 



•00 

47r sin^x • o^cilaj 





1 4 foo 

0"j, ( Vo) = K'l^Vo ^^^ TT sin^ 2x . a c?a. 







The definite integrals in the last two expressions are pure numbers depending on 
%2 alone. "^ We shall denote them by X^(%2) and \^^ {ux^ respectively. We thus 
havet 

so that, as in the case of molecules which are elastic spheres, Q'' is merely a constant 
multiple of 0' If we substitute the last expressions in the integrals for P, R, S (see 
(equations (42)) and remember that 



•00 



"y ^+1 ^ m-{- 



-, Von 



m 






we find that 



r 



(48) ^F„ = iK',,x'K)(^J 



/\3. 
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r 3^ 
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\ ^11-1 



iC^ —— g { O 
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^12 
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/I'll J- / 



yb, = -#,-(3-"-^^)(4 



2 



^2 — 2 5 



mo— I 



12 -^/ 



%2-l 



where 



K^, = (i^\'^i 



11 



?^n — 1 



^11 



* And on whether the molecules repel or attract, there being one value for each case, corresponding to 
each value of n. 

t These formulae are true whether the forces are repulsive or attractive. The only difference occurs in 
the value of the numerical constants X' and X". 
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Kii being the constant of attraction or repulsion between two molecules m in the 
same way as K^g is the constant for two unlike molecules. 

Thus h^ ki, h^ are in this case, as in the preceding case of § 13, numerical constants. 
The calculation of \^ and \^^ for various values of n from the definite integrals already 
indicated would afford an interesting theoretical investigation, but as actual molecules 
do not conform very closely to our hypothesis, the practical importance of the matter 
would hardly justify the labour. Maxwell calculated their values, however, by 
quadrature, in the special case ti = 5 (the forces being repulsive), and found that"^ 

(49) X'(5) = 2*6595, X''(5) = 1-3682. 



15. Rigid Elastic Spheres which Attract One Another. 

It is an undoubted fact that molecules attract one another at small distances (as is 
manifest from the force of cohesion, to take but a single example). By considering 
the effect of such forces in bringing about collisions between molecules (regarded as 
rigid elastic spheres) which would otherwise pass by one another, SuTHERLANDf had 
great success in explaining the variation * of viscosity with temperature. His 
treatment of the problem, while very suggestive, and forming an important 
contribution to the kinetic theory, laid no claim to rigour. 

In applying the present methods to the study of the same problem, we shall use the 
notation o-, o-^, 9S12, (}>i2 of §§ 13 and 14 ; ^^ir) will now represent the mutual potential 
due to the forces of cohesion. As before, we consider the path of the centre of the 
second molecule relative to that of the first. When the apsidal distance exceeds 
cr+cr^ no collision takes place, and the deflection 2x is given by the same equation as 
in § 14, viz., by 



(50) 



X = 



'no 





l-.^+^l.U^12(^)-^\2(-) 



V.^ 



^ 



-1/2 



dfj. 



When, however, the apsidal distance is less than o-+cr' a collision will take place, 
and the deflection 2x is twice the angle between the asymptote of the relative path 
and the radius vector from the centre of the first molecule to that of the second when 
the two are in contact {i, e, , for r = o- + 0-'). 

The differential equation of the relative path is 



2r^ 



1 fdr 



r^ \d0j p^Yq 



rMYo'+<t>'M-<i>'A^)}-i- 



The condition for a collision is evidently 



P 



' < ^-W^-{iVo^ + /x.(- + -) -^'i.(»)}- 



^ 



^ * Scientific Papers/ vol. ii., p. 42. X' (5) and X" (5) are, of course, tbe same as Maxwell^s 
Ai and A2. 

t ' Phil. Mag.,' 1893, (5), xxxvi., p. 507 
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Since the forces are attractive, — f^ is negative, so that 

dr 



(51) 



i>\A'T + <r')-<t>\A^) = -[H^ 



dr 



dr 



is positive ; we shall denote it by ^h^ for convenience. If we write p^ for 



cr-rcr 

^v7 



{v+Y,T, 



the condition for a collision becomes 



P <Py„- 



Next, by elementary geometry, the angle x, between the radii vectores r = cr + o-' 
and r = co ^ is given by 

, d9 

tan X — ''^ 



so that for p < py^ we have 



dr 



sin^X 



fl fdrV ,1-1 
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{rrwn^+Y,^) p,r 



sin^ 2x = 4 



^(^^ P' 
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1- 



Pv„ \ pv: 



Consequently we have 



^^2 C^o) = 4:7rVo sin^ X'Pdp = ^ttVo ^' -^ j9 (ip + 47rVo sin^ x-pdp, 

Jo Jo Py^ J« 

Vo 

where in the latter integral x is given by equation (50). We denote the latter 
integral by f\2 (^o)? the form of the function depending on the law of attraction 
between the molecules. 
Thus 



Q\. (Vo) - -Vo (--+-Oy;+Vo-) +4.V0 A, (Vo). 



^ 



Similarly we have 



Q"x.(Vo) = -kv, (^+'^01(^^+Vo-) +^Vo/^,(yo), 



where 



'^ 



.•X 



f\2 ( Vo) = I sin^ 2x . i) c?p. 



Substituting these values in the integrals P, R, S, we find that 
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where P\2 denotes the last integral on the preceding line. 
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With a similar notation for R'^^, &c., we find that 



Rt 






\hmm J I m-Vm ] 



^ ^ \hmml ) I ^ m+m^J 

The terms P\25 ^c., represent the effect of the attractive forces in deflecting the 
molecules without inducing collisions ; they are functions of the temperature in 

, and can be expanded m ascendmsf powers 01 >, a^.e., m aescenamff 

powers of 0. It may be shown that in each case they are of the form 



\(T-\-tT ) 



A3 / ^+m^ Y J A / hmm^ V _i_ "r / hmm^ V 4. 



\hnnmf / { \m-\-m! 



\m + m^, 



where 5 = 3, 4, or 5, according as we are considering P, B, or S. 
Hence the preceding formulae may be written 

2/m + m'Y/, , C12 . A 



F,,= .{.wr{ 



hmm! I 



1 + 



+ ^+•••1. 
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R' 



12 



•TT 






and so on, where we nave put (J10 = x :^ 7. 

SuTHEELAND, in his formula for the coefficient of viscosity, neglected all the terms 
of lower order in 6 than C12/0, and the success of the formula in representing the 
variations of viscosity with temperature (in the case of gases, though not in the case 
of vapours) seems to show that this procedure is legitimate under ordinary conditions 
in gases. This is equivalent to the neglect of the effect of the attractive forces in 
deflecting molecules without producing collisions. 

If we agree to neglect these terms our formulae become 
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(52) 
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* R denotes the universal gas-constant of the formula j?t; = (R/m) 

3 isr 2 
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where we have written 

n __ i_ ^^ 

In this case, therefore, the constants h are no longer independent of the 
temperature. It must not be forgotten that these last equations do not apply to 
vapours. 

Part III. — Discussion of the Theoretical Eesults, and Comparison with 

Experimental Data.'^ 

16. The principal formulae arrived at in the first part of this work are (33)~(4l); 
they are quite general and involve no hypothesis as to the mode of interaction 
between molecules (it being always understood that these are spherically symmetrical). 
Although these formulae indicate certain general laws, which require no knowledge 
of the nature of the encounters, for the most part the equations can be reduced to a 
useful form only when the quantities symbolically denoted by P\i, P^s? R\i? ^^? ^i? ai^d 
^2 are properly evaluated. This can be done only in certain special cases ; and in 
Part II. these calculations have been performed for three kinds of molecules, viz., 
rigid elastic spheres, point centres of force repelling or attracting one another 
according to the inverse n^^ power law, and rigid elastic spheres surrounded by fields 
of attractive force. The results of these calculations are given in equations (44)-(46), 
(48), and (52) respectively. 

When the formulae are thus completed, by comparison with experimental data we 
may find which of these three representations of the molecules best explains the facts, 
though, as all of them are somewhat artificial and ideal, it is not necessary that any 
one should excel the others in all respects. We shall find, however, that the third 
hypothesis (that the molecules are rigid elastic attracting spheres) is, on the whole, 
much the best of the three. It is, indeed, remarkable that so simple a mechanism can 
explain so much. 

As it might be considered that the choice of formula allowed by the general 
nature of our theory may possibly conceal errors due to some weakness in the 
foundations of the analysis, I shall first consider the cases wherein this latitude of 
choice is least, beginning with the case where there is none at all. I refer to 
equation (34), viz., 

which is a perfectly definite relation between three measurable quantities. I shall 

^ Part III. has been entirely rewritten in October-November, 1911; much new matter, not in the 
original Part III., has been added, this consisting chiefly (at Prof. L armor's suggestion) of further 
comparisons of the theory with experiment. Wherever possible I have used recent data, and I have in 
particular made great use of the results obtained by the pupils of Prof. Dorn, of Halle, from a well- 
directed series of experiments on viscosity, diffusion, and conduction. I am indebted to Prof. Dorn for 
the loan of many reprints and dissertations relating to these experiments. 
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next discuss the formulse for the conductivity and viscosity of mixed gases, as these 
only entail reference to the special law of interaction in the case of the numerical 
constants /c, k^^ and ^^ (whose values are, moreover, approximately known, as they do 
not differ very much with the special nature of the molecules). Lastly, I shall 
consider diffusion and the viscosity of a simple gas, as the corresponding formulae 
involve the most intimate molecular data, which can be measured only in very 
indirect ways. 

Any defect in analytical approximation involved in the theory is probably 
quite negligible under ordinary conditions of pressure and temperature ; the assump- 
tion made in § 4 concerning the law of distribution of velocities is the only place 
where error might arise, and both on theoretical and experimental grounds this 
appears to be inappreciable."^ Therefore, where no definite molecular hypothesis is 
involved, the agreement with experimental data should be perfect within the limits 
of experimental error. These conditions apply to equation (34) just quoted, and to a 
less extent to the formulae for conductivity and viscosity in mixed gases, where only 
slight weight attaches to the special law of force which must in their cases be 
assumed. 

Very different considerations apply to the formulae for diffusion and viscosity. 
These have great weight in indicating the best representation of the molecules by the 
manner in which they vary with the temperature, but if we use them to determine 
the molecular diameter which is involved in the expressions for D^2 ^^id yoc, we are 
treading on hazardous ground. Any discrepancy between the molecular diameter 
obtained from yoc and that obtained from D^^ must be ascribed rather to the artificial 
nature of the molecules postulated by the theory than to any defect in the theory 
itself From what modern electrical theories teach us, it appears extremely unlikely 
that any definite molecular diameter exists at all, even in the case of monatomic gases. 
Hence, while the formulae D^g and /x afford valuable independent means of determining 
the approximate dimensions of molecules, their interpretation must not be strained 
too far ; the agreement between the two sets of values may not necessarily lie within 
the limits of experimental error. 

While the present theory is strictly concerned only with monatomic gases, it is 
largely applicable to polyatomic gases ; for, as the molecules are in rapid motion, they 
must exert their actions equally in all directions during any short space of time, and 
will therefore behave very much as though they were spherically symmetrical. All 
general laws which hold good for monatomic gases may be expected to hold also for 
polyatomic gases ; and the results of experiment bear out this conclusion. Numerical 
agreement is not to be looked for, however, as the variable action of the polyatomic 
molecule will affect the numerical constants of the theory in taking the mean, just as 
would be the case if we had supposed all the molecules to be moving with the same 
speed, instead of different speeds varying widely about a mean value. In addition, 

■^ See the note on p. 483. 
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as the present theory does not deal in any way with internal molecular energy, it will 
not apply to the conduction of heat in polyatomic gases. 

17. The Coefficient of Thermal Conduction, 

Our expression for this quantity in terms ot the coefficient of viscosity and the 
specijfic heat at constant volume is 

which assumes that the molecules are spherically symmetrical, but is otherwise 
perfectly general. All previous writers on the kinetic theory have agreed on the 
conclusion that 

where f is a numerical factor ; but here agreement ends. Maxwell,'^ dealing with 
molecules which are point centres of force repelling one another according to the 
inverse fifth-power law, found as a result of his theory that /= f ; this is, of course, 
a very special case of the present theorem. Other writers, such as Clausius, 
Stefan, and 0. E. Meyer, have found values for / varying from 0'5 upwards, many 
of these calculations, however, being confessedly rough attempts, while Maxwell's 
theory attained a very high degree of accuracy. 

Maxwell's hypothesis is known not to be borne out by experimental facts (such 
as the variation of viscosity with temperature), and the theory of conduction which 
has hitherto found most acceptance is MeyepJs, which assumes that the molecules are 
rigid elastic spheres. Meyer's work was a valuable attempt at an exact treatment 
of the problem, but the method adopted (that used by Claijsiub and Maxwell in 
their early researches, but afterwards abandoned by Maxwell as being misleading) 
did not really allow of great accuracy. The expression for / which was arrived at 
involved a definite integral, which was calculated (using mechanical quadratures) by 
CoNRAU and Neugebauer.I" JeansJ improved the proof of Meyer's theorem, but 
his corrections did not affect the final result. The law obtained was 

& = 1*6027mC,; 

it has generally been considered that the difierence between this value of /, 
and Maxwell's value |, arose from the difierent nature of the molecules con- 
sidered. This view is in sharp contradiction to the theorem we have proved. 
Hence, since actual molecules do not conform to Maxwell's hypothesis, whatever 

■^ ' Scientific Papers,' vol. ii., p. 74. By a numerical slip lie gave the value of /'as |. The error was 
pointed out by Boltzmann and Poincar:^. 

t Meyer's ' Kinetic Theory of Gases' (English edition, 1899), Chapter IX. 
J Jeans' '' Dynamical Theory of Gases," Chapter XIII. 
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experimental evidence we can adduce in favour of the value /* = f will confirm that 
theorem. 

As both theories treat of monatomic gases, we must seek for evidence with regard 
to them. When this paper was written the only monatomic gas, of which I was 
aware, for which S- is known was mercury vapour, referred to by Meyer. ^ While 
a number of polyatomic gases obeyed his law fairly closely, this monatomic gas 
alone formed a striking exception. KooHf has determined yu for mercury vapour at 
273° C. and 380° C, and also at the much lower temperature 203° C. ; Schleier- 
macherJ has determined 9- at 203° C. ; these data, together with the theoretically 
calculated value of C^, lead to 3*15 as the value of f, Meyer raises some weighty 
objections against the accuracy of the data, viz., (i) that Koch's three values of /x 
show an unlikely amount of variation with the temperature ; (ii) that the conden- 
sation of saturated mercury in the capillary tube probably affected the determination 
of /x ; and (iii) that it is uncertain whether the vapour is completely monatomic at 
203° C. While the resulting value oi f is certainly unreliable, it is hard to conceive 
of the experimental errors being so great as to explain the difference between 
1*6027 and 3*15. So far as it goes, it tends rather to the support of the present 
theory. 

Quite recently, however, I have found that Sghwarze had disproved Meyer's 
theory nine years ago, by showing that /=2'5 for the true monatomic gases 
argon and helium. § As these determinations are important in this connection, a 
few details of them will be given. Schleiermacher's method, now generally 
accepted as the best, was used to determine 3-, heat being conducted by the gas from 
a heated platinum wire. The gases were very carefully prepared and purified,|| and 
in the reduction of the observations (which were made at two temperatures) due 
corrections were applied for the heat lost by conduction along the wire and for the 
temperature drop at the walls. The following values of S^o were obtained : — 

Argon 0*00003894, Helium 0*0003386, Air 0*00005690, 

the last named being determined as a check on the apparatus, and agreeing well with 

'^" ^ Kinetic Theory,' p. 295. 

t * WiED. Ann.,^ 1889, xxxvi., p. 346. 

I Ihid.^ 1883, xix., p. 857. 

§ W. SCHWARZE, 'Inaugural Dissertation,' Halle, 1902; 'Ann. d. Physik,' 11, p. 303, 1903; 'Phys. 
Zeitschrift,' 4, p. 229, 1903. 

Mehliss (' Halle Diss.,' 1902) made an earlier determination of ^ for argon by the Stefan- Winkelmann 
method, and found / = 2*44; this result was in such striking disagreement with Meyer's theory that 
Sghwarze investigated the matter very thoroughly, both for argon and for helium, by the more accurate 
method due to Schleiermagher. The result was quite unexpected, and neither observer seemed inclined 
to regard it as a confirmation of Maxwell's hypothesis, but rather as being due to a numerical defect in 
Meyer's theory — which view is nearer the truth, according to the present theory. 

II The helium contained \ per cent, impurity (probably of neon); the same mateiial was used by 
Sghultze ('Ann. d. Phys.,' 6, p. 303, 1901) to determine the value of \l used above. 
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Winkelmann's^ value 0*0000568 and MtiLLEn'st value 0*000056. MehlissJ obtained 
0*000038 for argon. 

SoHULTZE§ determined the following values of /uq for argon and helium, using the 
same materials as Schwarze : — • 

Argon 0*0002104, Helium 0*0001891. 

Using the following values of C^ : — • 

Argonir 0*0740, Helium|| 0*7142, 

Schwarze thus finds these values of/: — • 

Argon 2-501, Helium 2*507. 

I proceed to discuss briefly the case of polyatomic gases. It is immediately obvious 
that the present theory does not apply here, for C„ in actual fact differs widely from 
its value for the monatomic gas contemplated in the foregoing calculations (this is 
not to say that /is not equal to | for any polyatomic gas, but merely that our theory 
leaves the question perfectly open). The same remark applies to any theory which 
supposes the molecules to be devoid of internal energy — in particular to Meyer's 
theory. But the latter has derived all its support from polyatomic gases. Meyer's 
views on this point underwent some changes. In 1877"^^' he seems reluctantly to have 
accepted the theory (strongly upheld by Stefan and Boltzmann) that the internal 
and translational molecular energies travelled at different rates (the latter most 
rapidly), so that the conductivity would be less for a gas whose molecules possess 
much atomic energy than, for a similar gas with little atomic energy. 

In 1899,tt however, he held that the conductivity is the same for both kinds of 
energy, and supported this view by an unsound argument based on the law oi 
equipartition. This enabled him to assert that / is equal to 1*6027 for all gases, and 
so obtain all possible support from the data for polyatomic gases. More modern data 
would give much less support to the theory, as we shall see. But such disproof is 
unnecessary, for Schwarze's experiments conclusively show (i) that Meyer's value 

^ ' WiED. Ann.; 48, p. 180, 1893. 

t Ibid, 60, p. 82, 1897. 

J ^ Halle Diss.,' 1902. 

§ •' Halle Diss.,' 1901 ; ' Ann. d. Phys.,^ 5, p. 140, 1901, and 6, p. 302, 1901. 

II Calculated from the formula |R/ Jm ; m = 0*1439 for the helium used, and J was taken as 427. 
This leads to the above value of /, given in ' Ann. d. Phys.,' 11, p. 303, 1903. In ' Phys. Zeitschrift,' 8, 
p. 229, 1903, J was taken as 424, which made / = 2*490. 

^ Cv for argon is calculated from Cp (determined as 0*1233 by Dittenberger, 'Halle Diss.,' 1897) 
and 7 (determined as 1*667 by Niemeyer, cf. 'Smithsonian Physical Tables,' 1910, p. 232). Thus 
C^y - 0*0002104. 

^* 'Kinetic Theory of Gases,' 1st edition, 1877. 

ft Ibid , 2nd edition, English translation, 1899, pp. 291-296. 



CONSTITUTED OF SPHERICALLY SYMMETRICAL MOLECULES. 



465 



of / is theoretically unsound, and (ii) that the internal energy does not travel at the 
same rate as the translational energy. Hence the agreement of Meyer^s theory with 
experiment must be accidental. 

BoLTZMANN", who opposed Meyer's view, developed a theory of conduction"^' on the 
basis of Maxwell's hypothesis, taking into account the internal energy. He obtained 
the following relation between / and y, the ratio of the specific heats, 



/=¥(y-i). 

This reduces to Maxwell's law when y = f , but the formula, as we shall see, is 
not borne out by the experimental data. 

I have not attempted to work out a theory of conduction in polyatomic gases, and 
shall be content with pointing out how (in a general way, and with some marked 
exceptions, which may, however, be due to faulty data) / tends to be larger or 
smaller (while always less than |) according as a gas has less or more internal energy. 
The table below gives the values of \ for all the gases for which determinations 
are available, together with y, C^, /x, and the values of f calculated from them, and 
also / as calculated from Boltzmann's formula. The gases are arranged in increasing 
order of y, i.e.^ in diminishing order of ^, the ratio of internal energy to total energy. 



Gas. 


y.t 


Cy.J 


flQ X 10^. 


^0 X 10^. 


/ 
(observed). 


/ 
(Boltzmann). 


Ethane 


1*2191 


• 3024 


843 


4957 


1-94 


0-82 


Ethylene . . 






248 





•274 


944 


3958 




•53 





•93 


Carbon dioxide , 






300 





•1477 


1388 


3078 




•50 


1 


•12 


Nitrous oxide . 






304 





•1483 


1381 


3508 




•71 


1 


•14 


Methane. . . 






3162 





•45P 


104 


746^ 




•59 


1 


■19 


Ammonia . . 






3173 





•3976 


960 


4588 




20 


1 


•19 


Nitric oxide . 






394 





1665 


1680 


45P 




61 


1 


•48 


Oxygen .... 






402 





'1563 


1900 


57811 


^ 


95 


1 


•51 


Air ..... 






405 





•1695 


1721 


56910 




•95 


1 


•52 


Nitrogen . . 






405 





1738 


1670 


56911 




•96 


1 


•52 


Hydrogen . . . 






402 


2 


•427 


854 


387111 




87 


1- 


51 


Carbon monoxide . . 


1-409 


0-1730 


1628 


4998 


1-77 


1-53 



^ 'POGG. Ann.,' 157, 1876, pp. 457-469. The theory is partly empirical, being an adaptation of 
Maxwell's formula in which/ = |. Boltzmann states (p. 468) that the numerical coefficients would 
have to be altered if any other molecular hypothesis were adopted. Our theorem shows that this is not 
true, at any rate when y = f. 

t These values of y are taken from Jeans' ''Dynamical Theory of Gases " (pp. 220, 221), except where 
the contrary is indicated. 

1 LandolT and Bornstein's tables ; observed by Muller. 

2 'Smithsonian Physical Tables/ 1910; observed by Muller. 

8 Ibid. ; by WiJLLNER. 

I These values of C-y are taken from Jeans' treatise (p. 21 8)5 except where the contrary is indicated. 
They are due to Wiedemann and Wullner. 
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^ Calculated from y. 

^, ^ Calculated from C^ and y, as given in the Smithsonian tables. 

7 ZiEGLER, ^ Halle Diss./ 1904. Stefan ('Wiener Sitz./ 72, IL, p. 69, 1875) had found ^ox 10^-780 
for methane, and Winkelmann ('Pogg. Ann./ 156, p. 497, 1875) obtained the value 647. There are no 
previous determinations for ethane. 

s Determined by Winkelmann, as given in the Smithsonian tables. 

■ ^ Winkelmann (Meyer, p. 295). 

1^ ScHWARZE, loc. cit. ; Winkelmann's and Muller's determinations have already been mentioned. 
There are also experiments by Todd ('Roy. Soc. Proc.,' A, 1909, 83, p. 19) and Eokerlein ('Ann. d. 
Phys.,^ 3, p. 120, 1900), the latter being at low temperatures (for air, hydrogen, and carbon dioxide). 

^1 G-tJNTHER, 'Halle Diss.,' 1906. Also see Meyer's treatise for earlier determinations. 

The authorities for the above values of /xq will be given later, when we come to 
discuss the coefficient of viscosity. The table shows that the results agree neither 
with Meyer's nor Boltzmann's theory. The values of / for methane and ammonia 
stand out from the others^ which show a tendency to increase with y ; as all the 
values of / are less than f it would appear that the view of Stefan and Boltzmann 
is correct, that the atomic energy travels slower than the translational energy. 

In conclusion it may be remarked that the formula S- = /mO^, shows that S- will 
behave (as the temperature or pressure varies) in a way which can be predicted from 
the behaviour of /x and 0^ separately, if the equation is correct. Experiments have 
confirmed this^ and we shall therefore leave the discussion of these laws till we come 
to the simpler case of fi itself. 

18. The Coefficient of Conduction for a Mixed Gas. 

The expression we have obtained for this quantity enables us to determine it in 
terms of the coefficients of viscosity of the two pure gases forming the mixture, and of 
their coefficient of diffusion, provided that we know the law of interaction between 
the molecules. The latter is involved in the constants ^, k^^ and Aj^, but these do not 
vary very much with the law of force. In Part 11. of this paper we have determined 
their values in some special cases as follows : — ^- 

Rigid elastic spheres (46) — - 

/t —- 5 5 /i^i — "5 5 % 



'2 — 2 5? 



axwell's repelling molecules (48 and 49)— 

Jc =: 0771, h = Ip ^• 



— 7 
^^2 — 5? 



Attracting rigid spheres (52) — • 



7. - 3 1+39 Jl J. - 6 l + |Cl2 /g L _ 4 8 l+jOjd 



where (as we shall see in § 21) the constant C12 must be determined by experiments 
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on the variation of the coefficient of diffusion (of the gases concerned) with 
temperature. 

The formula for the coefficient of conduction in a mixed gas (equations 40 and 41) 
seems very compHcated, but the calculations from it in any particular case are quite 
simple. Fortunately we have experimental data for testing the law. Wachsmuth,'^ 
at Halle, has recently determined the conductivity of mixtures of argon and helium 
in various proportions ; the gases being monatomic, our formulae are properly 
applicable. Wachsmuth himself undertook the research, at Prof. Corn's suggestion, 
in order to determine f in the formula 9- = finQ^,^ taking the value of in from 
Tanzler's experiments (which will be discussed in the next section) on the viscosity 
of mixtures of argon and helium. He found that f so determined for the mixture 
was greater than f , approaching a maximum (about 4) when there was 60 per cent, 
of helium in the gas, and falling to f , as Schwarze showed, when either gas was 
eliminated. This fact is interesting, but, in the absence of any explanatory theory, 
does not lead to anything further. 

Wachsmuth also found that the observations could be represented to within 2 per 
cent, or 3 per cent, by the formula 

^- — '^ 4_ "^ 



(modelled on a similar formula for the viscosity of gases which we shall consider 
presently). As there were only four observations and two empirical constants the 
fact is not very remarkable ; when p = of course ^13 = 9-^, and when j9^ = we have 
^12 = ^. A much better agreement was obtained by allowing the empirical constants 
A and B (obtained by the method of least squares) to be imaginary. In this case the 
expression reduces to the quotient of a cubic by a quartic homogeneous expression in 
the variables j9, p^ ; as Wachsmuth remarks, the relation is good as an interpolation 
formula, but that is all. 

We proceed to determine the value of ^^2 from our formula. Taking oxygen as 
the standard gas (for which p^ = 0*001429 at 0° C and normal pressure) we have %v 
(for argon) equal to 1*224, and w^ — 0*125. The values of /xq ('^'•^•? M at 0° 0.) 
determined by Schultze {loc. cit. on p. 463), viz., 0*0002104 and 0*0001891 respec- 
tively, will be used ; also the previously given values of C^. The coefficient of 
diffusion D12 for argon and helium has been determined by Schmidt and LoNius 
(references given in § 21) ; reduced to 0° C. and normal pressure, in C.G.S. units it is 
0*650.t We now have all the necessary data for the calculation of 9i2 in the case of 
rigid spherical molecules and of Maxwellian molecules. For the case of attracting 
spherical molecules we need also the coefficient C12 ; this has not been determined, to 

^ 'Halle Diss.,' 1907; ' Phys. Zeitschrift/ 7, p. 235, 1908. The method was that due to Schleier- 
MACHER, and the apparatus was that previously used by SchwARZE and Gunther, aire ady quoted, 
t As explained in § 21, this value is uncertain to within 2 per cent, or 3 per cent, 

3 2 
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my knowledge, but the value 100 (intermediate between the values of 0^ for argon 
and helium) is near enough for our purpose, as the constant does not have much 
weight in the formula. 

From the equations (41) it is an easy matter to calculate the constants E, F, G, 
Ej, Fi, Gi ; we find the following values* for the three special kinds of molecules 
considered : — 



Rigid spheres- 



E = G-66, 
Ej = 3-93x10*. 



F = 2073, 
Fi = 376x10*, 



G = 10-60, 
Gi = 0-704 xlO^ 



Maxwellian molecules- 



E = 4-10, 
El = 2-42x10", 

Attracting spheres — ■ 

E = 5-38, 
El = 3-17x10*, 



F = 17-56, 
1= 2-65x10*, 



G= 10-65, 

0-706x10*. 



F = 19-0, 
Fi = 3-13x10*, 



G = 10-56, 
Gi = 0-700x10*. 



We have written Ei and Gi in place of l^ivf/m and Gw'/im' respectively. Our formula 
thus becomes 

^^12 ~ 2 \^vjy 



J2 



where (C^)i2 is given by equation (39a). The following table gives the values of p 
and p^ and the observed values of S-i2,t from Waghsmxjth's paper, together with S-^g 
calculated according to the above three hypotheses :- — ^ 

Thermal Conductivity of Mixtures of Argon and Helium. 



(argon). 


(helium). 


(Oz;)l2. 


^0 X 10^ 
(observed). 


^0 X 10*^ (calculated). 


Eigid 
spheres. 


Attracting 
spheres. 


Maxwellian 
molecules. 


1-000 

0-730 

0-546 

0-153 

0-0539 

0-0000 


0-0000 

0-270 

0-454 

0-847 
0-946 
1-000 


0-0740 
0-0984 
0-1257 
0-3094 
0-4952 
0-7142 


389 
741 
1077 
2320 
2939 
3386 


389 

675 

957 

2208 

2750 

3386 


389 
723 
1054 
2370 
2850 
3386 


389 
797 
1230 
2550 
2900 
3386 



^ These and nearly all the other calculations in Part III. have been performed with a slide rule, which 
is accurate enough for the purpose. 

t ^12 was determined at two tempera,tures ; we are, of course, using the value reduced to 0° C, as 
given by Wachsmuth, 
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Considering the nature of the calculation, the agreement is extremely good in all 
three cases (of course the values of 9- for the pure gases are taken from the formula 
34). The hypothesis that the molecules are rigid elastic spheres surrounded by fields 
of attractive force gives notably the best results, however, the agreement being 
within 3 per cent, in every case — which is little, if any, more than is introduced by 
uncertainty of the data. 

The above are the only available data for monatomic gases, so far as I know. 
Wassilewja^' has determined the conductivity of mixtures of hydrogen and oxygen, 
and, of course, we have the conductivity of air ; but as these are mixtures of 
polyatomic gases our formula does not apply to them, so that they will not be 
discussed here. 

In conclusion, it is hardly necessary to remark that our formula for the conductivity 
of a mixed gas reduces to that for a pure gas on putting p or p^ equal to zero. Since 
(as we shall see) D^s, the coefficient of diffusion, varies inversely as the total pressure, 
PqD^2 is independent of the pressure, as also are all the other quantities entering into 
E, F, G, Fi ; hence the coefficient of conduction of a mixture, like that of a pure gas, 
is independent of the total pressure. This has been experimentally verified in the 
case of air.t 



19. The Coefficient of Viscosity of Mixed Gases. 

The viscosity of mixed gases has been much studied, both theoretically and 
experimentally. It is especially interesting on account of the curious fact, first 
noticed by Graham ('Phil. Trans.,' 1846), that the addition of a moderate amount of 
light gas (hydrogen, in the case mentioned) to a much more viscous and heavy gas 
(carbon dioxide) may actually increase the viscosity of the latter. The same 
phenomenon was commented on by Maxwell in a Bakerian Lecture (1866), the 
gases receiving particular mention being air and hydrogen. 

The principal formulao which have been deduced for the viscosity of a mixture are 
due to Maxwell,J Puluj,§ Sutherland,!! and Thiesen.IF Of these only the first is 
based on an adequate proof, but as it is a particular case of my own formula (got by 
putting h = A2/A1 = 0771) it does not need separate discussion. Puluj's formula is 
based on an ingenious adaptation of the earlier theories of viscosity of a pure gas, and 
although, like those theories, it is only approximate, and in some ways not 

^ 'Phys. Zeitschrift,' 5, 1904, p. 737. 

t Meyer, * Kinetic Theory,' Chapter IX.; Todd, ' Eoy. Soc. Proc.,' A, 83, p. 19, 1909. 

J 'Scientific Papers,' vol. ii., p. 72. 

§ 'Carl's Eepertorium,' xv., p. 590; 'Wien. Sitzungsber.,' 1879, Ixxix., Abth. 2, pp. 97, 745; see also 
Meyer, ' Kinetic Theory,' p. 200. 
II 'Phil. Mag.,' 1895, xl., p. 421. 
U ' Verh. d. Deutsch. Phys. Ges.,' 4, p. 238, 1902. 
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satisfactory,"^ it nevertheless forms a useful interpolation formula. The expression is 

{mp^m'pj^' (Ptttttt^ 

Sutherland, by an argument which though not rigorous is very interesting and 
suggestive, arrives at the formula 

/x /x 



T—n- + -^ — ^—7—/ (Sutherland, Thiesen), 

which Thiesen also obtained by a different method of proof While, however, 
Thiesen left the formula as an empirical one, Sutherland strove to find expressions 
for A and B in terms of molecular constants. He found such expressions, but did so 
half empirically by a study of Graham's data. In the case of the gases which he 
considered, a very fair agreement was obtained. 
My own formula (39) can be written 

where E, F, G, F^ are given by equation (39) in terms of /x, jul\ D^g (the coefficient of 
diffusion), and a constant k which depends on the particular law of action between 
the molecules ; k is unknown, but may be expected to lie near or between the values 
already found for it in special cases (0'60 for elastic spheres, 0*771 for Maxwellian 
molecules). 

Of the above formulae, that by Puluj is the only one which is perfectly explicit. 
Thiesen's expression is completely empirical, and is useful only as an interpolation 
formula. PuLUj's relation shares the virtues and defects of the theory on which it is 
founded, and therefore prescribes a law of variation with temperature inconsistent 
with the facts for most gases. Thiesen's law, on the other hand, does not give any 
information concerning variation with temperature, and the constants A and B must 
be empirically determined for each case. 

The formula (34) or (53), as we shall see later, completely expresses the relation 
between jul^2 ^^d the temperature, but, being quite general, it specifies a difierent law for 
each molecular hypothesis ; and while, conversely, the determination of k and the law 
of temperature variation from experimental data may lead to further knowledge as to 
the best molecular hypothesis, this very generality gives the formula a semi-empirical 
character. 

On theoretical grounds it is desirable that the success of equation (34) as an 
interpolation formula should be tested. Excellent experimental material exists for 
the purpose. ScHMiTTf has lately given a resume of an extensive series of experi- 

^ Tlius it tacitly implies that the viscosity varies as the square root of the absolute temperature (since 
the theories on which it is based lead to this law). 

t ^ Ann. d. Phys.,' 30, p. 393, 1909. Full references to the original sources of the data are there given. 
ScHMiTT was apparently unaware of Sutheeland's formula, but, of course, it is the same as Thiesen's, 
if treated empirically. 
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ments on viscosity carried on at Halle, and has compared the results with the 

formulae of Puluj and Thiesen. The former agrees remarkably well with the 

observations, the residuals seldom exceeding 3 per cent. Thiesen's formula agrees 

even better still, the residuals rarely exceeding 1 per cent. This is not unnatural, as 

the constants A and B were found by the method of least squares to three or four 

significant figures. I have taken two of the longest tables given by Schmitt, and 

have roughly found the value of k (by trial and error) which makes the formula (34) or 

(53) agree most closely with the observations. The results are given in the following 

tables"^^ : — 

Viscosity of Mixtures of Argon and Helium at 15° C 



p 

(argon). 


(helium). 


/Xi2.10^ 

(observed). 


/X12.107 

(calculated) 
7^-0-692. 


1 • 0000 


0-0000 


2220 


2220 





•9507 





•0493 


2231 


2230 





9093 





■0907 


2243 


2240 





8571 





1429 


2253 


2249 





8074 





1926 


2266 


2259 





•7705 





■2295 


2264 


2266 





6846 





3154 


2266 


2283 





6119 





•3881 


2303 


2294 


0- 


5337 





•4663 


2299 


2297 





2915 





•7085 


2280 


2292 





1921 





•8079 


2226 


2246 


0-0000 


1 • 0000 


1966 


1966 



Viscosity of Mixtures of Oxygen and Hydrogen at 15° C. 



(oxygen). 


P 
(hydrogen). 


/X12. 10'^ 
(observed). 


i 

/X12.107 

(calculated) 
^-0-65. 


0-0000 


1-0000 


878 


878 


0-0521 


0-9479 


1092 


1076 


0-0878 


0-9122 


1191 


1191 


0-1561 


0-8439 


1359 


1370 


0-3333 


0-6667 


1674 


1672 


0-5678 


0-4322 


1877 


1876 


0-8126 


0-1874 


1992 


1975 


0-9555 


0-0445 


2014 


2007 


1 • 0000 


- 0000 


2014 


2014 



"^ In the first table the experimental constants used were as follows : /xxlO *^=^2220, /x'xlO^ = 1966, 
w = 1-224, id - 0-125 (relative to oxygen), po (at 15° C.) - 0-001353, D12 (at 15° C. and in C.G.S. 
units) = 0-705 (from Schmidt's experiments, cf. p. 478). The values of E, F, G, Fi were 7-897, 6*897, 
1^069, and 2-776 x 10* respectively. 

In the second table the values of fxxW were 2014 and 878, w - 1, w' ^ 0-0629, po = 0-001353, 
D12 (from Jackmann's experiments, c/. p. 478) — 0*760, at 15° C. and in C.G.S. units. The values of 
E, F, G, Fi were 11*34, 9*676, 1*049, and 5*205 x 10^ respectively. 
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The sum of the residuals between columns 3 and 4 is 77 in the first case, 54 in 
the second'^ ; more accurate determinations of h would probably reduce these 
figures, but the agreement is already within the limits of experimental error. It will 
be seen that the values of ^ which have been found, viz., 0*69 for argon and helium 
and 0*65 for oxygen and hydrogen, are intermediate between the values corresponding 
to the elastic-sphere theory and Maxwell's hypothesis, a result which confirms our 
theory. It is impossible to tell what special molecular structure these numbers 
indicate until the integrals of Part II. of this paper are worked out for other cases ; 
it may be that a very slight modification of one of the hypotheses there considered 
would explain the figures. 

The above tables, and others which might be given, show that our formula agrees 
well with observation, with a value of h accordant with the theory ; its success 
is therefore more significant than that of Thiesen's formula, with its two 
empirical constants. It may be noticed that Thiesen's expression can also be put 
into the form (53), with the speciality, however, that a relation exists between the 
three independent constants of (53) reducing them to two. Nevertheless, this 
agreement in functional form is sufficient to explain the success of Thiesen's 
expression as an interpolation formula over the limited range (O, l) of j) and p^. It 
may be concluded that, as far as the latter purpose goes, Puluj's formula, though 
not theoretically well founded, is the best (when the values of /x and /x' at the 
temperature considered are assigned), because it is very simple in form, quite definite, 
and sufficiently accurate for most purposes ; while Thiesen's formula requires the 
knowledge of two empirical constants, and my own formula is not so well suited for 
numerical calculation, f 

The study of the variation of the viscosity of mixed gases with the temperature is 
best deferred till the analogous question has been discussed for simple gases ; the 
subject will therefore be briefly discussed in § 22. 

20. The Coefficient of Viscosity of a Simple Gas. 

The formulae which have already been discussed do not contain any reference to 
the internal structure of the molecule, so far as this could be avoided, but are 

* The corresponding figures for Thiesen's formula are 55 and 51, and for Puluj's formula 126 and 
187. Of, K. SCHMITT, ^ Ann. d. Phys.,' 30, p. 393, 1909; pp. 408, 406. 

t The form of the expression (53) is well suited to explain the phenomenon mentioned at the beginning 
of this section. It shows that /xi2 is really a weighted mean between /x, /x' and F/Fi, the weight depending 
on the ratios E : F : CI and varying with p and p. For mixtures of CO2 and H we have F/Fi > /x > /x' ; 
also as E : G, or 1 + Jm/tv : 1 + Jm'/iv is large, on account of the largeness of w/w\ /x' is given very little 
weight in the mean except when p is very small. Therefore, for moderate values of the ratio p/p', the 
value of /X12 lies between F/Fi and /x, and may exceed the latter considerably. 

The same thing may be noted to a small extent in the above argon-helium table. From the figures 
already given, we have F/Fi = 2482 x 10""^, which is greater than 2220 x 10"^, the viscosity of argon at 
the temperature considered. 
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expressed as relations between directly measurable quantities. This enabled us to 
verify that our conclusions were borne out by the results of experiment. The 
equations might have been written otherwise, however, so that the formulae should 
express 3*, ^-ig? M12 '^^ terms of unknown molecular data, which might then be found 
from the experimental values of S- and ^12^* there would still be some check on the 
theory by the comparison of the values of the molecular constants found from 
different formulae. The reference to these constants was avoided by substitution of 
JUL and Di2 into the formulae, which were thus rendered much simpler ; no such 
simplification can be made in the expressions for jm and D^g, and as these involve the 
molecular data in the least complicated form, they are much the most suitable for 
the purpose of determininef the diameter of molecules. 
Our expression (33) for the viscosity of a simple gas is 



JUL 



^2w Y 



4h^mP \hfn/ ttR 



ff 






B^^u is a function of the temperature only, so also is h (since ih = 1/R0); and m, the mass 
of a molecule, is a constant for any gas. Hence we have obtained a perfectly general 
proof (for the case of a monatomic gas) that the coefiicient of viscosity is independent 
of the pressure and density of a gas, and depends only on the temperature. This 
remarkable law was first discovered by Maxwell. 

If we substitute the values of R^^n, which we have already found in special cases 
(see equations 45, 52, and 48), we get the following special laws : — 



(54) 



(55) 



— 5^ 1 1_ 



5m I R ^\^^^ 



Q4:^ir ' ^{2hm) ' (T^ 64(r^^7r\m 



JUL 



5m 



n ^Y' 1 



64(T^y^7rWV ^ , On 



(elastic spheres), 



(attracting spheres). 



1 + 



d 



(56) 



n+S 



^ =: Ae2(^-i) 



(centres of force oc l/r"^), 



where in the last equation we have 






m 



3 +■ 



23-4/(^-1) 



x''(^)r(4 



2 \ \4Knm, 



<mj 



n—1^ 



(the force between two molecules at distance r being 'Kum^jr^), and X^^ {n) is a number 
depending only on n^ having two values for each one of n^ according as the force is 
attractive or repulsive. 

The first equation shows that if the molecules are simple elastic spheres the 
viscosity varies as the square root of the absolute temperature (this law is due to 
Maxwell) ; the second shows that when the spheres also attract one another the 

VOL. OOXI. — ^A. 3 P 
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viscosity varies as 0^'^(l + C/0)~\ as was first discovered by Sutherland^ ; while the 
last equation shows that, if the molecules attract or repel one another according to 
the inverse n^^ power law, then /x is proportional to the (n + 3)/2 (7^— l)*^ power of Q. 
It is interesting to notice that Lord RayleighJ" predicted this law of variation from 
a consideration of dimensions alone ; the present formula is complete, with an exact 
analytical expression for the numerical constant. If we write /xq for the viscosity at 
0° C, the three laws just indicated may be written 

71 + 3 

/. _ (eY' M ^ (e\^' 1 + c/00 ix _ fe_Yn-i) 



Mo 



.Oj ' ^Q \6o/ l + C/0 /xo \9oJ 



When the molecules repel one another according to the inverse fifth power of the 
distance, the last equation becomes 

which was obtained by Maxwell. 

Experiment shows that the second formula (generally known as Sutherland's 
formula) agrees far better with the actual facts than the others. Maxwell's two 
hypothetical laws, /x oc 0^''^ and julozQ^ are obeyed (even approximately) by very few gases. 
The law /x oc 0' (where 6^ = (n + 3)/2 (^— l)) represents the variation much better in most 
cases, but even this applies only over a small range of temperature, after which a new 
value of 71 is required. On the other hand, Sutherland's formula applies universally 
so far as it has been tested.J It will be sufficient to mention the cases of helium and 
hydrogen. Schmitt§ gives tables showing that the law agrees with the observations 
of these gases from —60° C. to 185° C. ; below —60° C. the agreement ceases to be 
good. Thus this evidence tends to confirm the hypothesis that the molecules behave 
for our purpose like attracting spheres. 

The value of C has been shown in Part II. (pp. 457-460) to be a multiple of the 
mutual potential of two molecules in contact with one another. This was shown by 
Sutherland in his original paper, but as his mathematics was intentionally only 
approximate, his numerical constant is not correct. || The correction may be of 
importance when the theory of the constant C is further developed. IT It must be 
remembered always that 1 + G/9 is only the beginning of an infinite series of powers 

* ' Phil. Mag.; 1893 (5), xxxvL, p. 507. 

t " On the Viscosity of Argon as affected by Temperature," ' Scientific Papers,' vol. iv., p. 452, and 

p. 481. 

I We are referring to gases under normal conditions ; when the critical point is approached, the law 

ceases to be valid. 

§ 'Ann. d. Phys.,' 30, 1909, p. 399, where references to the original sources will be found, as well as 

data for other gases. 

II This does not affect his discussion (loc. cit) of the law of attraction between molecules, as he was there 
concerned only with the relative values of C for different gases. 

^ It is interesting to notice that Eankine ('Phil. Mag.,' January, 1911, p. 45, and *Eoy. Soc. Proc.,' 84 
p. 1 90) has found that for several gases C is proportional to the critical temperature. 
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of 1/0, and that later terms become important in the case of vapours, though for 
ordinary gases these, the first two, terms appear to be enough. 

As we have seen, the above laws of variation with temperature given by our theory 
have been announced before. The difference between our formulae and the earlier 
ones chiefly consists (as regards viscosity) in the numerical constants. We have 
already seen, in the case of the coefficient of conduction, to what grave errors 
(y= 1*60 instead of 2*50) the older elastic-sphere theory may lead. The expression 
given by Meyer for the viscosity differs considerably from the present one ; Jeans^ 
made an undoubted improvement in the older theory by allowing for the persistence 
of velocities after molecular collisions. The fact that this important consideration had 
been previously overlooked illustrates the imperfections of that theory. As it is, the 
correction for persistence of velocities is itself of uncertain amount t ; however, by its 
means the discordance between the formula of the older theory, and that given by 
the present theory for the case of elastic spherical molecules, is reduced almost to 
10 per cent. J 

From equations (54) and (55) we proceed to determine what evidence the various 
hypotheses lead to, as regards the size of molecules. Those formulae are strictly true 
only for monatomic gases, but can be applied to polyatomic gases with much more 
safety than could the formulae for the conductivity of gases, where the interna.1 
energy (which here plays little or no part) was very important. The formula (56) 
could also be used to determine the relative force constants of the molecules when 
assumed to be point-centres of force, but as this is of less interest we shall not trouble 
to do so. 

The said equations contain R/m, which can be determined from data as to the 
density of a gas at a given temperature and pressure. A table of values of R/m for 
several gases is given by Jeans (p. 113, loc, cit). To determine o-, the radius of the 
molecule, we also require to know m, or, since p = pm, we require to know i/, which is 
AvoGADRo's constant. Jeans used the value 4x10^^, the best one then available; 
but recent researches! agree in indicating 2*77 x 10"^^ as the correct values, which we 
shall therefore use here. In the following table are given the values of juq, C, and 
o-xlO^ for some gases, the latter being calculated according to the two hypotheses 
considered. It will be noticed that the radii of the molecules on the hypothesis that 
they attract one another is less than that when no such forces exist ; of course, in the 

* "Dynamical Theory of Gases,^' p. 238, p. 250. 

t Ibid.^ p. 250, line 13, where it is assumed for the sake of simplicity that the excess of momentum 
above that appropriate to the point of collision goes in equal proportions to each molecule. 
J In our notation, Jeans' formula (581) may be written 

•88 m /R^Y^^ 



which is '896 times our own expression. 

§ Sutherland ('Phil. Mag.,' 1909, February, p. 320) quotes Eutherford as the authority for this 
figure, and mentions that it is in good agreement with Planck's value 2*80 x 10^^. 

3 P 2 
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latter case the apparent size includes part of the extension of the field of force of the 
molecule. In selecting the values of juq and C I have preferred to use the mean 
of recent determinations (where these exist) rather than to go back to Graham's or 
Obermayeb's experiments. 



Gas. 


1 


C. 


o" X 


108. 












Elastic-sphere 


Attrac ting-sphere 








theory. 


theory. 


Air 


17211 


111-5 


1-830 


1 - 542 


Hydrosren .... 


8542 


75-6 


1-332 


1- 


180 


Oxygen ..... 


19003 


130-3 


1-781 


r 


469 


Nitrogen .... 


16724 


111-7 


1-842 


1- 


553 


Nitrous oxide . . . 


138P 


167 


2-267 


1" 


787 


Nitric oxide . . . 


16805 


195 


1-869 


1- 


429 


Carbon dioxide . . 


1388^ 


249 


2-259 


1" 


636 


Carbon monoxide . . 


16287 


156 


1-868 


1 


490 


Chlorine ..... 


12805 


199 


2-654 


2' 


019 


Ethylene .... 


9905 


249 


2-389 


1 


•730 


Helium 


18858 


75-3 


1-065 





•943 


Argon ..... 


21079 


162 


1-790 


1 


•508 


Neon ...... 


298110 


56 


1-266 


1 


•154 


Krypton ..... 


233410 


142 


1-622 


1 


•317 


Xenon 


210710 


252 


2-408 


1-738 



1 Authorities: The value of f^o is a mean of 1714 ('Landojlt and BOrnstein's Tables'), 1713 (Hogg, 
*Proc. Amer. Acad./ 40, 18, p. 611, 1905), 1733 (Breitenbaoh, *Ann. d. Phys.,' 5, 1901, p. 166), 
1730 (SOHULTZE, 'Ann. d. Phys.,' 5, 1901, p. 140), 1736 (Tanzler, 'Verh. Deutsch. Phys. Gesell.,' 8, 
p. 221, 1906), and 1702 (Grindley and Gibson, 'Roy. Soc. Proc.,' A, 80, p. 114); the yalue of C is a 
mean of 119*4 (Breitenbach), 113 (Sutherland, 'Phil Mag.,' February, 1909, p. 320), 111-3 (Ray- 
LEIGH, ' Scientific Papers,' vol. iv., p. 452, p. 481), and 102*5 (Grindley and Gibson). 

2 /xo is the mean of 864 ('Landolt's Tables'), 857*4 (Breitenbach, loc, dt.), 841 (Markowski, 'Ann. 
d. Phys.,' 14, p. 742, 1904); C is the mean of 71-7 (Breitenbach), 72*2 (Rayleigh, loc. dt.\ 83 (Mar- 
kowski, loc, cit), 79 and 72 (Sutherland, 'Phil. Mag.' (5), 36, p. 507, 1893, and February, 1909, p. 320). 

3 The mean of 1873 (Landolt) and 1926 (Markowski, loc, cii.) ; C from 128*2 (Rayleigh, 'Roy. Soc. 
Proc.,' 62, 1897; 66, 1900; 67, 1900), 138 (Markowski), 127 and 128 (Sutherland, loc. cit.). 

4 Graham and Obermayer {cf. Meyer's 'Kinetic Theory') 1660, 1670, and Markowski, 1675; 
C from 115, 113 (Markowski), 109, 110 (Sutherland, loc. cit). 

5 /xo from Jeans' treatise, p. 251 (due to Obermayer); C from Sutherland, loc. cit 

6 Breitenbach's value of /^o, and his and Sutherland's values of C (239, 258) ; both as previously 
cited. 

^ /xo from 'Landolt's Tables'; from. Sutherland, 'Phil Mag.,' February, 1909, calculated from 
' Landolt's Tables.' 

8 //.o from 1891 (Schultze, loc. cit), 1887 (Sghierloh, 'Halle Diss.,' 1908, cited by Schmitt, 'Ann. d. 
Phys.,' 30, 1909, p. 393), 1879 (Rankine, 'Phil. Mag.,' January, 1911, p. 45); and C from 80*3 
(Schultze), 72*2 (Rayleigh), 78*2 (Sghierloh), 70 (Rankine), 76 (Sutherland, loc. cit). 

^ p.Q from 2104 (Schultze, loc. cit), 2114 (Sghierloh, loc. cit), 2102 (Rankine, loc. cit); from 169*9 
174-6, 142 (by the preceding three authorities), and 160 (Sutherland, loc. cit). 

10 Rankine (loc. cit). 
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As we shall see, independent values of the radii of molecules may be obtained from 
the coefficient of diffusion, which we proceed to consider. 

21. The Ooefficient of Diffusion. 

As the coefficient of diffusion of a gas into itself can be expressed much more 
simply than that of one gas into another, we will discuss it first, though it cannot be 
measured experimentally. 

The expression obtained in Part I, was 

Q /^\l/2 1 



1^11 =-~o (^^) 



A2/ ' ^ F 






Since v is proportional to the pressure, it follows that D^ varies inversely as the 
pressure, the other factors being functions of the temperature alone. If we substitute 
for P^n from equations (45), (48), and (52), and eliminate the molecular constants by 
means of the expressions found in the preceding section for the corresponding values 
of ^, we get the following results : — ■ 

Dii = I ~ (rigid elastic spheres), 

P 

Dn = I : w \ ( 3 — - — 7 ) - (centres of force KmV""*), 
K {n) \ n—lj p ^ ^' 



1 -T"j: 



2C 
3 " 



Du = % 7T ~ (attracting spheres), 

d 

where in the last case is the constant in Sttthbblakd's law of viscosity. All these 
formula agree in showing that D^ is a numerical multiple of ju/p, the factor generally 
being a function of the temperature, though with the first two hypotheses above, the 
factor is a constant. When, in the case of the second hypothesis, we take n = 5 (as 
Maxwell did, the force being repulsive), the formula becomes 

Dji = 1'543 ~ (Maxwell's fifth-power law), 

p 

using the values already given for Ai and A^. The extension to general values of n 

is interesting, if only for the unobvious character of the factor (3 ~), 

\ n—lj 

The best value of the numerical factor hitherto obtained, on the hypothesis that 

the molecules are elastic spheres, is due to Jeans (see his treatise, p. 273). After 

allowing for persistence of velocities, he deduces that 1*34 is its approximate value; 

this is in very fair agreement with f, the result of the present theory. 
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The expression obtained in Part I. for the coefficient of diffusion between two 
gases is 

Since 

1 m ^ne 

where p^ is the total pressure, and since v and v do not occur elsewhere in the 
expression, it follows that, according to the present theory, the coefficient of diffusion 
is independent of the relative proportions of the gases in the mixture. On this point 
Stefan and Maxwell were at issue with Meyer,'^ whose theory predicted a large 
change in the coefficient of diffusion as the proportions of the gases were varied. 
Several experiments t have recently been made at Halle to test the rival theories. 
The pairs of gases taken were 0~H, 0-N, H-N, A-He, COg-H. The ratio of the 
components was varied from 1 : 3 to 3 : 1 (roughly), and the results showed systematic 
differences in the values of D12 which reached 8 per cent, in extreme cases. Though 
part of this variation may be due to experimental errors, it is undoubtedly the case 
that appreciable variations of D12 exist, or, at any rate, variations in the values of D12 
as derived from the accepted theory of the experiments. This theory is of rather an 
elaborate nature, and takes no account of the boundary conditions ; any attempt at a 
revision of the theory of diffusion should include an examination of the theory of the 
experiments by which the diffusion is measured ; and it is by no means impossible 
that faults in the accepted theory of the experiment may account for the small 
variations of 5 per cent, or so in the determinations of D^s. At the same time it must 
be remembered that the value of D12 here arrived at theoretically is based on the 
assumption that the square of the velocity of diffusion, together with small deviations 
from Maxwell's law of distribution, are negligible. If these had been taken account 
of (and it is easy to do so) part or the whole of the outstanding discrepancy would 
probably be explained. As the variations are so small, however, I have not troubled 
to carry out these calculations. 

The experiments emphatically disprove Meyer's theory, which predicts changes of 
20 per cent, or more in Di2j and there can be no doubt that the general principles of 
the theory as laid down by Stefan and Maxwell are much nearer the truth than 
those advocated by Meyer. A modification of Meyer's theory by Gross,J which 
reduced the amount of the variations predicted by Meyer, has also been shown to be 
incorrect, since the actual variations of D^g, though similar in amount to the theoretical 
values deduced by Gross, are in the opposite direction, 

* See his ' Kinetic Theory,' Chapter VIIL 

t R. Schmidt, 'Halle Diss.,' 1904, and 'Ann. d. Phys.,' 14, 1904, p. 801; E. Deutsch, 'Halle Diss.,' 
1907; Jackmann, 'Halle Diss.,' 1906; and LONIUS, 'Halle Diss.,' 1909. All these results are summarized 
by LoNlus, ' Ann. d. Phys.,' 29, 1909, p. 664. 

I G. Gross, ' AVied. Ann.,' 40, p. 424, 1890 
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The special forms of the coefficient of diffusion, corresponding to the three molecular 
hypotheses considered in Part II., are as follows : — ■ 

-n - _3 A /2 1 /I , 1\" f^Y 1 B> f^T ^^igid elastic\ 

JJ12-1F V ^(^+^')2^^ + ^'; [J ,^ pKeJ \ spheres > 



«+3 



'' 8 K'A'(n)rf3 2__\lU w'j'p,dj p ' O,' v, \of force/' 

-n - ^ A /I 1 A , lY'" (B^T 1 ^ f^Y' 1 /attracting 

1^12 - 16 V ^(^+a-OA^ ie^V \ J ^a P W l + Gje \ spheres 

where after substitution from equations (45), (48), and (52) we have put the formulae 
into a more convenient form as above ; w and w^ are, as before, the specific gravities 
relative to a standard gas whose density at pressure po and temperature Oq is p^ ; vq is 
AvoGADRo's constant, at the same pressure and temperature. It may be noticed 
that Maxwell,^ by applying the methods of his fifth-power-law theory to the case 
of a gas whose molecules were elastic spheres, obtained an expression | times as great 
as (but in other respects identical with) the first of the equations just given ; as he 
did not give any details, it is impossible to say how his error arose — probably by a 
numerical slip. 

The following equations simply express the law of variation of Dig with temperature, 
as shown by the above three equations : — 

Di2/(Di2)o = {O/OqY'^ (rigid elastic spheres), 



71 + 3 ^1 



^12/(^12)0 = (d/Oo)^^''-^^ (point-centres of force), 

Di2/(Di,)a = {e/e,f" { 1 -I- Gje,}/{ l -l- G,,/6} (attracting spheres). 

So far as I am aware, the variation of the coefficient of diffusion with the 
temperature has never been properly examined experimentally ; the values of D12 are 
generally found only at one temperature, or two at most, and this is insufficient to 
decide between the second and third of the laws just given. From the analogous 
behaviour of viscosity, as affected by temperature, however, the third law is probably 
most nearly true, and the values of C12 have been worked out by Sutherland on this 
hypothesis.'!* In view of the importance of the constant C12, it is desirable that 
further experiments be made on the variation of the coefficient of diffusion with 
temperature. 

"^ ' Nature,' vol. viii. ; ' Scientific Papers,' vol. ii., p. 343. 

t The present theory shows that C12, the temperature constant for diffusion, is a different multiple of 
the mutual potential of two molecules in contact from that in the case of Cn, the temperature constant for 
viscosity. This fact was not indicated in Sutherland's investigation; see * Phil. Mag.' (5), 38, p. 1, 1894, 
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As in the case of the coefficient of viscosity^ we proceed to give the values of the 
molecular radii (according to the two theories of elastic spherical molecules^ with and 
without attraction) calculated from the coefficient of diffusion. From the latter we 
get (T+(r\ the sum of the radii of molecules of each kind considered. By taking three 
pairs of the same three gases, like O— H, N-H, H— N, we can derive the actual values 
of a- in each case, which may then be compared with the values obtained from the 
coefficient of viscosity. I have thought it preferable, however, to compare the values 
of (T-\-a obtained from diffusion with the sums cr-\-a of the values obtained from 
viscosity (on the corresponding molecular hypothesis). The values of the coefficient 
of diffusion (which are in C.G.S. units, and for the temperature 0"" C.) are taken from 
the Smithsonian tables (1910 edition) ; they are uncertain (partly for reasons already 
discussed, partly for other reasons) to within about 5 per cent. ; different authorities 
give values differing by from 5 per cent, to 10 per cent. ; I have taken the mean, 
where two values are given in the table cited. The constants G12 are taken from a 
table given by Suthebland'^ ; as they depend only on two values of D^^, and as C 
occurs in such a form as to be largely affected by experimental errors, these values 
are not well determined, though they are probably the best available. In several 
cases I have not been able to find the value of C. 





■ Coefficient 


Tempera- 
ture 
coefficient. 




(cr + cr'; 


) X 108, 




_ Elastic 


-sphere 


At tracting-sph ere 


Gases undergoing diffusion. 


of 


theory. 


theory. 




diffusion. 
















^^12' 


From 


From 


From 


J^'rom 








diffusion. 
2-99 


viscosity. 


diffusion. 


viscosity. 


Air-hydrogen 


0*661 


74-5 


3-16 


2-65 


2-62 


Air-oxygen ........ 


0-1775 




3-42 


3-61 






Carbon dioxide-air ...... 


0-138 


250 


3-73 


4-09 


2-70 


3-18 


Carbon dioxide-carbon monoxide . 


0-136 


— 


3-79 


^ 4 


■13 






Carbon dioxide-nitrous oxide . , 


0-0983 


380 


4-19 


4" 


53 


2-71 


3-42 


Carbon dioxide-oxygen . . . . 


0-1802 


187 


3-23 


4" 


04 


2-49 


3-10 


Carbon monoxide-ethylene . . . 


0-101 




4-61 


4" 


26 






Carbon monoxide-hydrogen . . . 


0-642 


70 


3-03 


3' 


20 


2-70 


2-57 


Carbon monoxide-oxygen . . . 


0-183 


113 


3-38 


3 


'63 


2-84 


2-96 


Hydrogen-carbon dioxide . . . 


0-538 


115 


3-29 


3 


'59 


2-76 


2-72 


Hydrogen-ethylene 


0-486 


— — ' 


3-48 


3 


72 






Hydrogen-nitrous oxide .... 


0-535 




3-30 


3-60 






Hydrogen-oxygen ...... 


0-679 


100 


2-94 


3*11 


2-52 


2-55 


Nitrogen-oxygen ...... 


0-174 


117 


3-47 


3-62 


2-87 


3-02 


Argon-helium T ....... 


0-650 




2-53 


2-85 







Considering the possibility of occasional exceptionally large errors due to the 

^ *PhiL Mag.,' 1895, 40, p. 421; 1894, 38, p. 1. 
t R. Schmidt, 'Ann. d. Phys.,' 14, 1904, p. 801. 
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number of experimental data entering into each of the above determinations of 
molecular radii, the agreement shown is very good. It is notably better on the 
hypothesis that the molecules attract one another than on the ordinary elastic sphere 
theory (the large discrepancy in the case of CO^-N^O is probably due to the value of 
Ci2, which is larger than the value of 0^ for either gas, and this is somewhat improbable). 
In the case of the latter theory the value of a-k-a-' determined directly from diffusion 
is in almost every case less than that obtained from, separate observations of viscosity ; 
there is a slight tendency for the same thing to show itself with the values deduced 
from the other theory. This is probably to be attributed to the artificial nature of 
our hypotheses ; the effective distance of reversal of motion between two different 
molecules may not really be half the sum of the effective distances for two pairs of 
like molecules, though the rigid sphere hypotheses assume this. The third hypothesis 
we have considered does not do so ; but a comparison of diffusion and viscosity on 
this hypothesis would only give us the relation between the force constants between 
the molecules of one kind and different molecules, while the interesting question is the 
one we have just tested, viz., how nearly alike will a molecule behave towards another 
of the same kind, and towards one of a different kind. 

The following table^^ may be. of interest, showing how the values of cr obtained 
from Van deb Waals' law agree with those obtained from viscosity (on the 
attracting sphere theory) : — 



Gas. 



(Tan deb Waals' 
theory). 



Hydrogen . . 
Nitrogen . . 
xxir .... 
Carbon dioxide 
Helium . . . 
Argon . . . 
Krypton . , 
Xenon . 



(TX 10^ 

(attracting sphere 
theory). 



1-15 


1- 




76 


1- 




•64 


1- 




•70 


]• 




•16 


0- 




43 


1- 




56 


1- 




71 


1- 



18 
55 
64 
49 
94 
51 
32 
74 



22. The Variation of the Viscosity of Mixed Gases with Temperature. 

This question was deferred till after the similar questions regarding the viscosity 
of simple gases, and diffusion, had been discussed. The formula we obtained for (j.^^, 



^ 



The values of o- for the first four gases are calculated from values of h (in Van der Waals' equation) 
given by Jeans (/Dynamical Theory of Gases/ p. 141) from the experiments of Van der Waals and 
EosE Innes; but the value of v has been taken as 2*77 x 10^'^ instead of 4*0 x 10^^ as used by Jeans. 
The values of h for the last four gases are taken from a paper by RuDORF, 'Phil. Mag./ June, 1909, 
p. 795. 



482 MR. S. CHAPMAN ON THE KINETIC THEOEY OF A GAS 



was 






where 



G: ^ ^ 



i^ + 'l^VoI>i2" 



We have abeady seen that, on the elastic sphere hypothesis, and the point centre 
of force (KigmmV"'') hypothesis, 

where 5 = -J in the former case, and varies with n in the latter case. Hence 



P0D32 ^ 



s 



so that G, 1/m, and l//>c^ all vary with the temperature in the same way (we are 
supposing that the value of n is the same for both kinds of molecules). Moreover, we 
have seen in Part II. that Iz is independent of the temperature in these hypotheses. 
Hence the numerator of the above fraction is independent of 0, while the denominator 
varies as Q~\ In other words, the coefficient of viscosity of the mixed gas behaves in 
exactly the same way as that of the component gases, as the temperature varies. 

The case is rather different on the attracting sphere hypothesis ; remembering the 
expressions already found for k and /x and D12 it may easily be seen that 



D4-E0+F0=^ 



where for any given mixture A, B, C, D, E, F are independent of the temperature. 
Over a limited range of it is possible to find a constant K such that 



B + C0 A' 



D+E0+F0^ K + 

is extremely small, so that we may write 

Ai0^'^ 



Ml2 



l + K/0 



showing that according to the attracting-sphere theory Sutherland's formula holds 
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for mixtures as well as for simple gases — indeed, the formula has been more often 
tested for air than for any simple gas. Sohmitt"^ and others have shown that 
experimental results fuUy bear out this conclusion, t 

^ * Ann. d. Phys./ 30, p. 393, 1909; Tanzler Q Yerh. Deutsch. Phys. Ges.,' 8, 12, p. 122, 1906). 

t Thiesen (' Yerh. DeutscL Phys. Ges.,' 8, 12, p. 236, 1906) states that Sutherland's formula and 
his own law (see § 18 of this paper) are incompatible, although Tanzler's experiments demonstrate the 
contrary. 

Note added February 15, 1912. 

Later investigation shows that the neglect of the terms of higher degree than the third, in the 
expression for the function F (p. 441), is not easily justifiable by rigorous analysis. The omitted terms 
have been considered by Enskog ('Phys. Zeitschrift,' xii., 58, 1911); his work suffices to show that the 
formulae in this paper are correct except as regards a factor which in general is a function of the 
temperature (only), and in particular cases is a numerical constant. The latter cases include the case of 
elastic spherical molecules (as pointed out by Enskog), and also of point centres of force varying inversely 
as the n^^ power of the distance. Thus Lord Eayleigh's theorem (p. 474) can be rigorously established. 

I have endeavoured to form a numerical estimate of the correction factors to be applied to the formulse, 
but the analytical difficulties involved have thus far proved unsurmountable. The expressions obtained 
do not converge rapidly, and the calculations are very laborious ; but the problem is now perfectly definite, 
and improved analysis may remove the difficulties. 

The present theory must therefore be regarded as approximate only. In conjunction with Enskog's 
work, however, it provides rigorous proofs of the relations connecting viscosity and conductivity with the 
pressure and molecular diameter and mass, and in special cases also with temperature; the relations 
themselves are well known in the case of simple gases, but not in that of compound gases. The numerical 
constants are not rigorously determined, and are subject to correction ; but the agreement with experiment 
seems to show that the approximation is a good one. 

Lastly, it should be remarked that the theory of diffusion is unaffected by the terms omitted, so that 
the above statements do not apply to it. 



